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1.0 The lidar model 


1.1 Description 

The coherent lidar model was developed to allow the assessment of various space-based coherent 
Doppler lidar instruments under consideration at NASA Marshall Space Right Center as part of 
the Laser Atmospheric Wind Sounder (LAWS) program. This section describes the theory used in 
the model. 

1.1.1 Principle of operation 

A Doppler lidar wind sounder measures the wind velocity at a remote location. It does this by 
transmitting an optical pulse which has a narrow frequency spread. This optical pulse is backscat- 
tered from aerosol targets (for the optical frequencies under consideration here) in the atmo- 
sphere. These aerosol particles are moving with the velocity of the local wind field and therefore 
impart a Doppler shift to the backscattered light. The backscattered light is collected at the lidar 
receiver and mixed with an optical local oscillator to produce a heterodyne beat signal on a detec- 
tor. This heterodyne frequency (less any offset between the LO and transmitter laser) is propor- 
tional to the velocity difference between the lidar and the aerosol target. For a stationary lidar this 
difference is due solely to the aerosol velocity along the line of sight of the lidar. For a moving 
lidar, the lidar velocity along the line of sight of the instrument must be removed to obtain the 
aerosol velocity. 

1.2 Lidar geometry 

Figure (1.1) shows the geometry for a coherent lidar at an altitude, h above a spherical earth of 
radius, R*. The line of sight of the instrument makes an angle, 0 O , called the nadir angle with the 
vertical from the earth to the lidar. The nadir angle at the surface is On,, and is given by:- 

0ns = asin (ir) (U) 

' e ' 

where R1 is the length of the perpendicular from the line of sight to the earth’s center and is given 

by:- 

R1 = Hsin(0 ) where H = R +h (1.2) 

n c 

The range from the lidar to the surface is given by:- 

, 0.5 

R = (H +R;-2HR e cos(0 ns -0 n )) (1.3) 

For a target at an arbitrary altitude, z above the earth’s surface, the nadir angle at the target and 
range from the lidar are given by substituting (R* + z) for Re in equations (1.1) and (1.3) respec- 
tively. i.e 
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(1.4) 


e„(z) 



2 2 

R(z) = (H +(R e + z) -2H(R e + z)cos(0 nz -0 n )) (1.5) 



Figure (1.1) The geometry for a space-based lidar. 

1.3 The coherent lidar signal to noise ratio equation 

For an aerosol target at range, R with uniform backscatter coefficient, {3 over the dimensions of 
the transmitted beam, detector noise dominated by Local Oscillator (LO) shot noise the Coherent 
Laser Radar (CLR) SNR is given by: 

SNR - Trar 1 ! 0 ,, 

8h\)BR 2 


where U L is the laser energy, T t is the intensity transmittance of the transmit optics due to apertur- 
ing and obscuration, T tar is the intensity transmittance of the transmit optics due to losses 
(absorption, reflection) that are independent of the beam cross-section position, T A is the one-way 
atmospheric attenuation, D is the effective diameter of the transmit telescope, T^r is the intensity 
transmittance of the receive optics, c is the speed of light in a homogeneous atmosphere, h is 
Plancks constant, v is the laser frequency, B is the signal bandwidth of the detector, 'n mix (some- 


, 110 n mi , 

1 + r 

2T„P„ . 
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times referred to as r| h ) is the CLR mixing efficiency assuming perfect alignment and perfect 
polarisation matching, T R is the intensity transmittance of the transmit optics due to aperturing 
and obscuration , p 0 is the transverse coherence length of the backscattered optical field at the 
CLR, and T| is given by: 


^ ^mix^A^P^HQE 


(1.7) 


where r| A is the change in mixing efficiency due to misalignment between the back propagated 
local oscillator (BPLO) and the signal beam, T) P is the BPLO/signal beam polarisation mismatch 
factor and T| H q E is the detector heterodyne quantum efficiency. The transverse coherence length of 
the backscattered optical field is given by: 


P 0 (R) 



( 1 . 8 ) 


where H=2.914383 , k=2nl"k (X is the optical wavelength) and C„ 2 (x) is the refractive index struc- 
ture constant at range x from the CLR. 


1.3.1 Transmitted pulse energy 


The term U l T t T tar represents the energy transmitted by the lidar system. The loss mechanism 
due to the optics has been split into the two terms T t and T tar to allow various lidar geometries to 
be considered. T tar is simply the product of the intensity transmission efficiencies of the individ- 
ual optical components in the transmit path. T t is more complex and its value is dependent on the 
geometry of the lidar. Various authors [1-4] have shown that different receiver geometries and 
aperture truncations of the lidar transmit/receive beams will vary the heterodyne efficiency of the 
lidar system. Rye and Frehlich [3] analysed the optimum aperture truncation required for three 
lidar geometries. A summary of their results for Gaussian beams and large detector sizes (no 
detector aperturing effects) is presented in Table (1.1). 


Table (1.1) Optimum Lidar Geometry Parameters 


Geometry 

Transmitter 

Receiver 


System- 

Antenna 

Efficiency 


*"trunc 

T t 

** trunc 

t r 

mmmm 


Wang [2] 

1.763 

0.955 

1.763 

0.955 

0.42 

0.401 

Rye 1 ( 1 1 

1.626 

0.929 

1.626 

0.929 

0.457 

0.424 

Rye 2[ 1 1 

1.736 

0.951 

1.192 

0.758 

0.461 

0.438 


For a given physical diameter, D the 1/e 2 spot size, c is related to the truncation radius, r, n>fW . by 


o 


D 


r, Jl 
trunc 


(1.9) 
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and T t is given by:- 


T x = 1 - e 


-£)' 


( 1 - 10 ) 


Local Oscillator 


Secondary 

Optics 



Geometry I: An untruncated local osci- 
llator is combined with the receiver sig- 
nal. This is the geometry proposed by 
Wang and most commonly used in a 
practical lidar system. 


Photodetector 


Beamsplitter 


Local Oscillator 


Geometry II: The local oscillator is 
truncated in secondary optics desig- 
ned to optimise the BPLO field by 
avoiding its truncation at the rec- 
iever aperture. Designated Rye 1 in 
Table 1.1. 


Receiver 

Aperture 


Apertured 

Secondary 

Optics 



Photodetecton 


Beamsplitter 


Local Oscillator 


Receiver 

Aperture 



Photodetector 


Geometry III: An alternative to II. 
This is the theoretically optimum 
geometry. The local oscillator is trun- 
cated directly by the receiver aperture 
but the BPLO is again untruncated. 
Difficult to implement in a practical 
scheme as the local oscillator would 
require a beam diameter greater than 
that of the collecting telescope. Rye 2 
in Table 1.1. 


Figure (1.2) The three receiver geometries dicussed by Rye and Frehlich[3]. 
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The system-antenna efficiency, T) s is given by:- 





( 1 . 11 ) 


where A R and Ap are receiver and transmitter aperture areas respectively. The lidars under consid- 
eration here (monostatic lidars) have a common transmit/receive aperture and therefore A R =A r 
and the system-antenna efficiency is simply:- 

= ^mix T T 0- 12 > 


1.3.2 Atmospheric attenuation 

Atmospheric attenuation of the lidar beam arises from molecular and aerosol absorption and scat- 
tering. The Air Force Geophysical Laboratory has an atmospheric transmission/attenuation model 
based on the Hitran molecular database. This model or an equivalent can be used to generate 
atmospheric extinction/km as a function of altitude for each of the wavelengths of interest. The 
lidar model developed here uses these extinction as a function of height profiles to calculate the 
extinction between the lidar and the target aerosol. The model divides the atmosphere into 0. 1 km 
thick layers from 0 km to 20 km. Above 20 km the extinction is insignificant and can be ignorred 
for the purposes of this model. The atmospheric extinction, T A from the lidar to a target in the 
atmosphere is then given by:- 


T a = e " (1.13) 

where n is the number of atmospheric layers from 20 km down to the target, a„ and 1 0 are the 
extinction coefficent and path length respectively for layer n. The path length, 1 0 in a layer 
between altitudes zl and z2 is simply:- 

l n = R(z2) -R(zl) (1.14) 

where R(z) is given by equation (1.5). 


1.3.3 The refractive index structure function. 


The WPL-37 refractive index profile [5] is an experimental refractive index structure function 
profile which has been used in lidar modeling [6]. It can be approximated using: 


C n (z) = 


0 

3.66 x 10 


-17 


1.59 x 10 12 x z (m) _4/3 


7.36 x 10 


-14 


for 


z > 20km 
20km >z> 3km 
3km > z > 10m 
10m > z 


respectively (1.15) 
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Fairall and Frisch of NOAA have developed an analytic model of the refractive index structure 
function which takes into account location, season and local conditions!?]. This is an extensive 
model and the reader is referred to their technical memorandum for details, however the following 
typographical corrections were mentioned by the authors and should be noted: 

Equations 19 and 24c on pages 7 and 8 of reference [7] should read: 

L = T/ (Kg) (ujT.) 


and: 

h t ( 5 ) = 10(1 -^)- 3 / 2 ( 2-^) _1 

respectively and Table 1 of reference [7] should be altered and added to as follows:- 




Colorado 



Florida 



Winter 

Summer 

Winter 

Summer 



Night 


Night 

Day 

Night 

Day 

Night 

Po 

1 

10 

2 

10 

0.35 

10 

0.25 

10 

a 

0.4 

0.9 

0.2 

0.9 

0.2 

0.9 

0.2 

0.9 

u« 

0.5 

0.1 

0.35 



0.1 

0.35 

0.1 

R„„ (Win 2 ) 

300 

-100 



350 

-100 

700 

-100 

Zj(m) 

1200 

100 

2000 

100 

2000 

100 

2000 

100 

z, (m) 

7000 

7000 

12000 

12000 

9000 

9000 

14000 

14000 

T 0 (K) 

278 

263 

298 

288 

288 

283 

303 

293 


Table (2.2) Parameters for the Fairall and Frisch C„ model. 


Additionally the refractive index structure function for the atmosphere of the ocean trade winds 
(e.g. off Florida) can be calculated using the Florida summer day parameters with a, changed to 

0.93. 

Fortran code to generate the experimental WPL-37 c n ^ profile and the Fairall and Frisch models 
was developed (See “Appendix A - Refractive Index Structure Function Calculation” on 
page 23.). Appendix A also includes plots of these profiles. It can be seen that the model profiles 
and the WPL-37 profile are in general agreement. For the purposes of the lidar analyses under 
consideration here, variations in the refractive index structure function have small effects on the 
lidar SNR and so it was decided that the model would use the WPL-37 profile as a generic profile 
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1.4 The efficiency term 

This consists of the four terms r| A , % and T| H q E . The mixing efficiency (assuming perfect 
alignment), is geometry dependant and is given in Table (1.1). The variation of the mixing 
efficiency as a function of the misalignment angle [9] between the BPLO and signal beams is 
shown in Figure (1.3), where D is the telescope diameter and A. is the optical wavelength. 



Figure (1.3) Mixing efficiency as a function of JtD/X. 

As has been indicated previously the maximum mixing efficiency is dependant on the geometry 
used and so for the purposes of the lidar model the alignment efficiency, T| a was obtained by nor- 
malising the data used to plot Figure (1.3) by dividing by the maximum mixing efficiency. The 
mixing efficiency for a system with misalignment is then given by taking the mixing efficiency 
due to geometry considerations (Table (1.1)) and multiplying by tj a . 

The values of atmospheric backscatter, (5 typically used are based on experimental field results 
using an assortment of lidars with assorted polarisation states. It may therefore be assumed that 
losses due to polarisation are explicitely accounted for in the use of experimentally derived |J val- 
ues and therefore additional losses due to polarisation are likely to be small. However for com- 
pleteness, the polarisation efficiency, % , although set to unity, is included in the model 
formulation. 

The detector heterodyne quantum efficiency, T| H q E is the efficiency with which the detector con- 
verts coherent photons into photoelectrons at the maximum bandwidth of the signal. For a detec- 
tor on a spacecraft orbiting at -7 km/s this corresponds to a bandwidth of ~ 1.4 GHz at 9 pm and 
-12 GHz at 2 pm. The heterodyne quantum efficiency is usually significantly smaller than the DC 
quantum efficiency. 
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1.5 The velocity estimator 

The model uses a parameterised maximum liklihood estimator! 10] provided by Dr. Rod Frehlich. 
The model provides values of the standard deviation of good estimates, g and the fraction of bad 
estimates, b as functions of <I>, the effective number of photo-electrons detected for fixed values of 
M, the number of samples and ft, the effective number of independent samples per estimate. The 
values of these parameters are given in Table (3.3). The variation of g and b as functions of d> are 


Parameter 

Complex 

Real 

Description 

T 

L o 



Observation time per estimate. 

V 

T max 



Maximum velocity seen by processing algorithm. 1 

CO 



Spectral width of return signal. 

SNR 



Signal-to-noise ratio at processing bandwidth. 

T 

*s 

W.v„„) 

V(8.V„„) 

Time between samples. 

M 


t/t s 

Number of samples/estimate. 

ft 

CdM.T s 

qiM.T s 

Effective number of independent samples/estimate. 

0> 

SNR.M 

SNR.M/2 

Effective number of photo-electrons per estimate. 

Table < 

'3.3) Paramei 

ters used by the velocity estimator. 


a. It should be noted that there has been some confusion over the definition of V m „ , for a 
processing window of -V to +V (ie maximum wind speed of V but with direction away 
from or towards the lidar), V,,^ has been defined as 2.V or as V by various people. In 
this document V,,^ is always defined as V. 


given by the empirical relations: 


-(*) = 
to 


1 + 


<t> 

So 




-6 


and : 


(1.16) 


b(<I>) 



(1.17) 


Values of x< go* £. 8, p, b 0 , a and y for various values of M and ft are listed in Appendix C - Lidar 
Performance Calculation. 
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1.6 The lidar scan pattern 

The model considers three different scan patterns, a simple conical scan where the whole tele- 
scope is pointed off nadir and rotated (the original LAWS scan technique), a conical scan using a 
rotating nadir pointing telescope with a front refracting element to scan the beam and a line scan 
mechanism similar to the French BEST scan pattern. 

The shot pattern on the ground of the two conically scanned systems is the same, however the per- 
formance is different due to a reduction in the effective aperture of the telescope for the refracting 
element scanner. The effective telescope diameter, D c ff is: 

D eff = Dcos(0 n ) (1.18) 

and it is this value that is used for the diameter in the SNR equation( 1 .6). 



Figure (1.4) Geometry of effective beam diameter dependence on nadir angle. 

The line scan technique is discussed in Appendix B- AEOLUS Telescope Geometry on page 32. 

1.7 Results of the model 

NASA Marshall Space Flight Center is in the process of conducting studies on a small satellite 
coherent Doppler wind lidar known as AEOLUS (Autonomous Earth Observing Lidar Utility 
Sensor). The model was used to assess the performance of four concepts under consideration for 
this instrument. Two of the concepts use a 9 pm C0 2 laser device whilst the remaining two con- 
cepts use a 2 pm Ho:Tm:YLF laser device. For each of the two wavelength options, two different 
shot scanning schemes were considered - a dual-telescope bi-perspective line scan scheme and a 
conical scan using a rotating refractive wedge element. The laser device characteristics were cho- 
sen as those most likely to be available for a space mission without incuring a heavy development 
cost penalty. The parameters used by the lidar model for each of these systems is given in Table 
(4.4). 
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Parameter 


Design 1 | Design 2 j Design 3 Design 4 Units 


Orbit 



Orbit height 


Inclination angle 


Laser 


Wavelength 

9.11 

9.11 

2.06 

2.06 

pm 

Pulse length 

0.5 

0.5 

0.5 

0.5 

ps 

Pulse energy 

400 

400 

200 

200 

mJ 

PRF 

20 

20 

10 

10 

Hz 


Orbit duty cycle 


Optics 


Telescope diameter 


Nadir angle 




'Hp 

— 

BPLO/signal beam align- 
ment tolerance 


Receiver 


Type 

— 

Geometry 


Scan 

— 

Type 


Telescope rotation rate 


System margin 

— 

Target 


0.5 

0.5 

0.5 

0.5 

30 

30 

30 

30 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 . 

1 

1 

1 

1 

13.91 

16.07 

3.15 

3.64 



Complex Complex Complex Complex 



Dual-look Wedge Dual-look Wedge 



Table (4.4) AEOLUS designs analysed. 
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Parameter 

Design 1 

Design 2 

Design 3 

Design 4 

Units 

Altitude 

300 

300 

300 

300 

m 

Max. horizontal velocity 

30 

30 

30 

30 

m/s 

‘ 

Vertical range resolution 

1000 

1000 

1000 

1000 

m 


Table (4.4) AEOLUS designs analysed. 


The BPLO/signal alignment angle tolerance was chosen to give a value of 0.5 for T| a . This enables 
the tolerance requirement for the two wavelengths to be compared. An additional parameter listed 
in the table is the system margin. Traditionally the predicted and measured performance of coher- 
ent Doppler lidars have differed. The system margin is a 3 dB allowance to ensure that the design 
of the instrument will be adequate for the task. The atmospheric model used for the atmospheric 
extinction coefficient, T A was a mid-latitude summer clear aerosol model. This is the only atmo- 
sphere model that is included in the model at the moment and the data for this model were sup- 
plied by Coherent Technologies Inc., of Boulder, CO. Section shows some generic screen shots of 
the model (not for one of the four cases listed above). The performance of the lidar is generally 
characterised by two parameters, the backscatter value at which b=0.5 (1.17), i.e the velocity 
algorithm correctly determines the correct velocity 50 % of the time and by g (1.16) the standard 
deviation of the good velocity estimates about the true velocity. For the systems under consider- 
ation here the value of g at b=0.5 is always adequate for the task and so the performance of these 
systems can be characterised by the single parameter, (3(50 %) the value of aerosol backscatter at 
which b=0.5. 



Design 1 

Design 2 

DesignS 

Design 4 


(3(50 %) 

7.8 x 10' 9 

1.1 x 10' 8 

5.1 x 10' 8 

6.8 x 10* 8 

/(m-sr) 


Table (5.5) 

It should be pointed out that the AEOLUS design is still being refined. 

1.8 Summary 

A model of a space-based coherent lidar has been developed and is being used to assess the per- 
formance of various instrument designs. 
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1.9 Screen shots of the model 


The following pages show screen shots of various views of the model together with samples of the 
input parameter dialogs. 


ur 
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One benefit of having the input menus accessible from the button bar is the ability to change the 
lidar parameters from anywhere within the model. Changing the laser wavelength in this example 
would cause the extinction plot to automatically be recalculated and displayed for the user when 
the OK box on the input dialog is selected. 
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Samples of the other input dialog boxes. 



Msnm. page is 
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2.0 Appendix A - Refractive Index Structure Function Calculation 

2.1 Program structure 

Program CN2CALG - main routine, allows user to select c n 2 profile desired, determines routines 
to be used and calls them, saves results to output file called c n 2 .dat. Output is in three columns, 
first column is altitude in km, second column is c, 2 and third column is c n 2 . 

Subroutine CLRSCRN(I) - is a utility subroutine that scrolls the output screen by I lines. Called 
by the main program. 

Function WPLCN(Z) - is a function for calculating the WPL-37 c„ 2 profile. For this c n 2 case, the 
c, 2 profile is not generated. Called by the main program. 

Function CT2(Z,HP,P0,PBLZ,RNET,USTAR,ALPHA,BETA0) - is a function for calculating the 
c, value for a given altitude for all cases except the WPL-37 profile. This value is then used to 
calculate the value of c n 2 (except for WPL-37). Called by the main program. 

Function CN2(Z,HP,P0,CT,T0,ZT) - is a function for calculating the c n 2 value at a given altitude 
for all cases except the WPL-37 profile. Called by the main program. 

Function CTPBL(Z,PBLZ,ALPHA,RNET,BETAO,USTAR,T,DAY) - is a function for calculating 
the c, 2 value in the planetary boundary layer. Called by function CT2. 

2.2 Program code 
2.2.1 Program CN2CALG 

PROGRAM CN2CALG 

C 

IMPLICIT NONE 

C 

REAL *4 HP, PO, PBLZ , RNET , USTAR , ALPHA , OLDCT , 

♦BETAO , TO , ZT, ETA, MAXZ , DZ , Z , CN2 , CT, CN, CT2 , 

+BETA ( 12 ) , ALPH ( 12 ) , U ( 12 ) , NETR ( 12 ) , ZPBL ( 12 ) , 

+TZ ( 12 ) , TEMPO ( 12 ) , WPLCN 
INTEGER I , N, DAY, ILOCAL, ISEAS, ITIME, INDEX, 

+WPL 

c 

PARAMETER (HP=6500 . , P0=1000 . ) 

MAXZ=2 0000 . 0 
N=500 
WPL = 0 
OLDCT=0 . 0 
C 

DATA BETA / 2 . , 10 . , 1 . , 10 . , 0 . 25 , 10 . # 0 . 35 , 10 . , 0 . 25 , 0 . 25 , 0 . 25 , 0 . 25/ 

DATA ALPH / 0 . 2 , 0 . 9 , 0 . 4 , 0 . 9 , 0 . 2 , 0 . 9 , 0 . 2 , 0 . 9 , 0 . 93 , 0 . 93 , 0 . 93 , 0 . 93 / 

DATA U /0.35,0.1,0.5,0.1,0.35,0.1,0.5,0.1,0.35,0.35,0.35,0.3s/ 

DATA NETR / 6 0 0 . , - 1 0 0 , 3 0 0 . , - 100 . , 7 00 . , - 1 00 . , 3 5 0 . , - 100 . , 7 00 . , 

>700 . , 700 . , 700 . / 

DATA ZPBL /2000 100 1200 100 ., 2000 ., 100 ., 2000 ., 100 2000 . , 

>2000 . ,2000 . , 2000 . / 

DATA TZ /12000. , 12000. ,7000. ,7000. , 14000. , 14000. ,9000. ,9000. , 
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ono nno non non 


c 

c 


c 


+ 14000. , 14000. , 14000. , 14000 . / 

DATA TEMPO /298 . , 288 . , 278 . , 263 . , 303 . # 293 . , 288 . , 283 . , 303 . , 303 
+ 303 . , 303 ./ 


CALL CLRSCRN ( 40 ) 


WRITE ( * , 
WRITE ( * , 
WRITE ( * , 
WRITE { * , 
WRITE ( * , 
WRITE (*, 
WRITE ( * , 
WRITE ( * , 
WRITE ( * , 
WRITE { * f 
WRITE ( * , 
WRITE ( *, 
WRITE ( * , 
WRITE ( * , 
WRITE ( * , 
WRITE { * , 
WRITE ( * , 


*)'This program calculates the refractive index' 

*) 'structure constant, Cn2 for various locations, ' 
*) 'seasons and conditions. The calculation is' 

*) 'based on NOAA Technical Memorandum ERL WPL-195' 
* ) "Diurnal and Annual Variations in Mean' 

*)' Profiles of Cn2"' 

*)' C.W. Fairall and A.S. Frisch' 

*)' Wave Propagation Laboratory' 

*)' Boulder, Colorado' 

*) ' March 1991' 

*) ' 1 

*) ' This program developed by:' 

* ) ' Gary Spiers ' 

*) ' Center for Applied Optics' 

*>' University of Alabama in Huntsville' 

*)' Huntsville, AL 35899.' 

*) 


10 CALL CLRSCRN ( 2 ) 

WRITE {*,*)' Please select one of these locations:-' 
WRITE ( * , * ) 

WRITE {* , * ) ' 1) Colorado' 

WRITE { * , *) ' 2) Florida' 

WRITE(*,*)' 3) Trade Wind Oceans' 

WRITE ( * , * ) ' 4) WPL-37 Profile' 

READ < * , * ) ILOCAL 


TRAP INVALID LOCATION 

IF { (ILOCAL.LT. 1) .OR. ( ILOCAL. GT. 4) ) THEN 
WRITE (*,* ) 'Please choose a valid location!' 
GOTO 1 0 

TRAP WPL EXPERIMENTAL PROFILE 


ELSEIF ( ILOCAL. EQ. 4) THEN 
WPL= 1 


TRAP INCOMPLETE OCEAN TRADE WIND PROFILE 
ELSEIF ( ILOCAL. EQ. 3) THEN 

WRITE (*,*)' At present this program only calculates a 
+'daytime Trade Wind Profile.' 

ISEAS = 1 
ITIME= 1 


DETERMINE SEASON AND TIME OF DAY FOR REMAINING CASES 
ELSE 

11 WRITE (*,*) 

WRITE (*, # ) Please select one of these seasons:-' 
WRITE (*,*)' 1) Summer' 

WRITE (*,*)' 2) Winter' 

READ ( * , * ) ISEAS 

IF ( ( ISEAS. LT.l) .OR. { ISEAS . GT . 2 ) ) THEN 
WRITE (*,*)' Please choose a valid season option!' 
GOTO 11 
ENDIF 

12 WRITE { * , * ) 
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WRITE (*,*)' Please select the time of day:-' 

WRITE (*,*)' 1) Daytime' 

WRITE ( * , * ) ' 2) Nightime' 

READ ( * , * ) ITIME 

IF ( ( ITIME . LT . 1 ) . OR . { ITIME . GT . 2 ) ) THEN 
WRITE (*,*) ' Please choose a valid time of day!' 
GOTO 12 
ENDIF 
ENDIF 

WRITE ( * , *) 'Calculating. . . ' 

IF (WPL.EQ.O) THEN 
INDEX = ( I LOCAL- 1 ) *4+ (ISEAS-1) *2+ITIME 
DAY=2- ITIME 
BETAO=BETA ( INDEX) 

ALPHA=ALPH ( INDEX) 

USTAR=U( INDEX) 

RNET=NETR< INDEX) 

PBLZ=ZPBL< INDEX) 

ZT=TZ ( INDEX) 

T0=TEMP0 (INDEX) 

END IF 

OPEN ( UNIT=9 , FILE= ' cn2 . dat ' ) 

DZ=ALOG10 ( MAXZ ) /N 
DO 15- 1 = 1 , N 
Z=10 . ** ( I*DZ) 

IF (WPL.EQ.l) THEN 
CT=0 

CN=WPLCN ( Z ) 

ELSE 

ETA= Z/PBLZ 

IF ((ETA .GE. 0.95) .AND. (ETA. LT . 1.0)) THEN 
CT=OLDCT 
ELSE 

CT=CT2 { Z , HP , P0 , PBLZ , RNET , USTAR # ALPHA f BETAO , 

@ TO , DAY , ZT ) 

ENDIF 

OLDCT=CT 

CN=CN2 ( Z , HP , P0 , CT , TO , ZT ) . 

ENDIF 

WRITE (9, 20) Z/1000,CT,CN 
15 CONTINUE 
20 FORMAT (3E15 . 9) 

STOP 

END 


2.2.2 Subroutine CLRSCRN 

SUBROUTINE CLRSCRN(I) 
INTEGER J , I 
DO 100 J=1 ,1,1 
WRITE { # , * ) 

100 CONTINUE 
RETURN 
END 
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2.2.3 Function WPLCN 


REALM FUNCTION WPLCN (Z) 

C 

IMPLICIT NONE 
REALM Z 

C 

IF (Z.LT.10.) THEN 
WPLCN=7 .36*10 . ** (-14 . ) 

ELSEIF ( (Z.GE.10. ) .AND. ( Z . LT . 3 0 00 . ) } THEN 
WPLCN= 1 . 59*10. **(- 12. )*Z**(-4./3.) 

ELSEIF ((Z.GE.3000.) .AND. ( Z . LE . 2 00 00 . ) ) THEN 
WPLCN=3 .66*10. **(-17.) 

ELSE 
WPLCN=0 . 

END IF 

RETURN 

END 


2.2.4 Function CT2 

REAL* 4 FUNCTION CT2 ( Z , HP , P0 , PBLZ , RNET , USTAR, ALPHA , 
@BETA0 , TO , DAY , ZT) 

C 

IMPLICIT NONE 
C 

REAL *4 HP , P0 , DTROP , CTT , CTS , PBLZ , RNET, USTAR, ALPHA, 
*BETA0 , TO , ZT, Z, P, T, CT, CTPBL 
INTEGER DAY 
C 

PARAMETER (DTROP=2000 . ,CTT=8 . E-5 , CTS=1 . E-3 ) 

P=P0 *EXP { -Z/HP) 

IF (Z.LT.ZT) THEN 
T=T0- (T0-216 . ) * Z / ZT 
ELSE 
T=2 16 . 0 
END IF 

IF ( ( Z .LT. (ZT-DTROP)) .AND. (Z . GE . PBLZ)) THEN 
CT2 =CTT 

ELSEIF ( ( Z .GE. (ZT-DTROP)) .AND. (Z .LE. ( ZT+ 

* DTROP))) THEN 

CT= (ALOGIO (CTS) -ALOGIO (CTT) ) * (Z-ZT) / (DTROP* 2) 

* + (ALOGIO (CTT) +ALOG1 0 (CTS) ) /2 
CT2 = 10 * * (CT) 

ELSE 
CT2 =CTS 
END IF 

IF (Z .LT. PBLZ) THEN 

CT2 =CTPBL { Z , PBLZ, ALPHA, RNET, BETAO , USTAR, T, DAY) 
ENDIF 
RETURN 
END 
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2.2.5 Function CN2 


REAL*4 FUNCTION CN2 ( Z , HP, PO , CT, TO , ZT) 
IMPLICIT NONE 

REAL* 4 HP, PO , Z , CT, TO , ZT, P, T 

P=P0*EXP ( -Z/HP) 

IF (Z.LT.ZT) THEN 
T=T0- (TO-216 . ) *Z/ZT 
ELSE 
T=2 16 . 0 
ENDIF 

CN2= ( (7 . 9E-5*P/T**2 ) **2) *CT 
RETURN 
END 


2.2.6 Function CTPBL 

REAL*4 FUNCTION CTPBL ( Z , PBLZ , ALPHA, RNET, BETAO , USTAR, T, DAY) 
IMPLICIT NONE 

REAL *4 Z , PBLZ , ETA, HSO , ALPHA, RNET, BETAO , TSTAR, RHOCP, USTAR 
* , L , KARMAN , G , WSTAR , T , THETA , HB , HTB , HT , RT 
INTEGER DAY 

PARAMETER (RHOCP=1100, G=9.81, KARMAN=0 . 4) 

ETA=Z/ PBLZ 

HSO= ( 1 . -ALPHA) *RNET/ ( 1 . +1 . /BETAO) 

TSTAR= -HSO/ ( RHOCP*USTAR) 

L= (T*USTAR* *2 ) / ( KARMAN*G* TSTAR) 

IF (DAY . EQ. 1) THEN 
WSTAR= ( ( -G/T) *USTAR*TSTAR*PBLZ ) ** ( 1 . /3 . ) 

THETA= -USTAR* TSTAR /WSTAR 

HB= 1 0 . * ( -PBLZ/L) ** (2 . /3 . ) * (ETA* (l.-7.*Z/L))**(-2./3.) 
HTB=8 .*(1.-0. 8* ETA) 

IF (ETA .GT. 0.95) THEN 
HT=10 . * (0.05) ** (— 3./.2) /1.05 
ELSE 

HT= 1 0 . * (1 . -ETA) **(-3./2.)/(2. -ETA) 

ENDIF 

RT=-0 . 2 * (1.+3 .2* (-L/PBLZ) ) 

CTPBL=THETA* * 2 * PBLZ* * (-2-/3.)* ( HB+RT*HTB+RT* *2*HT) 

ELSE 

CTPBL=TSTAR* *2 *Z* * (-2./3.)*(l.+2.7*(Z/L)**{2./3.)) 

ENDIF 

RETURN 

END 


2.3 Plots of output results. 
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3.0 Appendix B- AEOLUS Telescope Geometry 

This appendix was originally a memo released to the AEOLUS team on 9/29/93. 
In this document, nadir always refers to the spacecraft nadir. 

AEOLUS 


Telescope 



Figure (1.5) 3-D geometry. 

3.1 Angle definitions 

< DEB = 9 n , the angle between the line of sight and nadir (nadir angle). 

< AED = 0 R , the angle between the left/right component of the line of sight and nadir. 

< DEC = Sp, the angle between the forward/aft component of the line of sight and nadir. 
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Consider triangle DEB: 


E 



L.sin(On) B 


Figure (1.6) Triangle formed by line of sight and nadir. 

Thus we have DE and DB in terms of L, the slant range to the ground and 0 n , the nadir angle. We 
can fold the triangles ADE and CDE into the ground plane ABCD: 



Figure (1.7) Folding of triangles ADE and CDE into the ground plane ABCD. 
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From Figure (1.6) we obtain expressions for ED and DB in Figure (1.7) and from Page (1) we 
have < DEC = 0 F and < AED = 0 R . Simple trigonometry gives expressions for AD and CD. ABCD 
is a rectangle and so AB = CD. AD, AB and BD are related by Pythagoras’ theorem: 

AD 2 + AB 2 = BD 2 (1.19) 

substituting gives: 


L 2 .cos 2 (0 n ).tan 2 (0 R ) + L 2 .cos 2 (0 n ).tan 2 (0 F ) = L 2 .sin 2 (0 n ) ( 1 .20) 

dividing by L 2 .cos 2 (0 n ) gives: 


tan 2 (0 R ) + tan 2 (0 F ) = tan 2 (0„) 


( 1 . 21 ) 


Thus we have obtained an expression relating 0 n to 0 R and 0 P 


3.2 Optimum geometry for wind velocity measurement. 

Ideally we would like to resolve the wind velocity components along and perpendicular to the 
ground track with equal resolution. In Figure (1.7) these components are represented in the 
ground plane by DC and DA respectively and for them to have equal weighting we should there- 
fore arrange the geometry such that: 

DA = DC (1.22) 

Substituting values and simplifying gives: 

0r = 0f (1-23) 

Substituting back into Equation (1.21) gives: 

2.tan 2 (0 F ) = tan 2 (0 n ) (1 .24) 

rearranging to give 0 F as a function of 0 n : 

tan(0 F ) = tan(0 n )/ (1.25) 


3.3 Angle between forward and aft view opticalaxes 

For the dual telescope option, the angle, a between the forward and aft view optical axes is the 
angle between the forward and aft lines of sight. As the geometry is symmetric about the line per- 
pendicular to the ground track, this angle is given by 2.< AEB.Figure (1.8) shows the geometry 
for this calculation where B' is the aft-view point equivalent to B. We see that : 

sin(a/2) = cos(0„).tan(0 F ) ( 1 .26) 

substituting for tan(0 F ) from Equation (1.25) and simplifying gives: 

a= 2.sin' , (sin(0„)/V2) (1-27) 

For AEOLUS we have 0 n = 30 deg. and therefore a = 41.41 deg. and 0 F = 0 R = 22.21 deg. 
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Figure (1.8) Geometry for telescope axis separation and path length on surface of earth. 


3.4 Path length between forward and aft points on the ground. 

Figure ( 1 .8) also shows the center of the earth, C e and the arc of the earth’s surface B'A'B between 
the forward and aft line of sight ground points. We see that B'AB represents the chord linking 
these two points. From triangle ABC e we can determine the half-angle subtended at the earth’s 
center, y/2 as: 

sin(y/2) = L.cos(0 n ).tan(0 F )/R < ( 1 .28) 

The angle subtended by the earth’s arc is y and the length of the arc is given by R^.y, where y is in 
radians. This gives: 

track length = IRe.sin-UL.cosfOJ.tanfOpVRe) ( 1 .29) 

For the AEOLUS conditions of 0 n = 0.5236 (30 deg), 0 F = 0.3876 (22.21 deg) and Re = 6371315 
m we have a line of sight length. L of 407916 m for a 350 km satellite altitude, assuming a spher- 
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ical earth. This gives a track length of 288498m between the forward and aft line of sight ground 
points. The ground track velocity is 7303 ms' 1 and therefore the time to cover this distance will be 

39.5 s. This means that to overlap the forward and aft views on the same patch of ground, the 
instrument should look in the forward direction for 39.5 s before changing to the aft view. At a 
constant PRF of 10 Hz the seperation between pulses on the ground will therefore be 730.5 m, for 
a 20 Hz pulse rate this separation will reduce to 365.3 m. 

3.5 Effect of varying 0„ . 

The following plots show the variation of 0 R , 0 F and 2. a as the nadir angle, 0 n varies between 0 
deg. and 90 deg. 

The first plot shows the dependence of 0 R on 0 F for various 0 n . 



Figure (1.9) The dependence of 0 R on 0 F for various nadir angles, 0 n . 

All of the remaining figures use the optimum geometry, 0 R = 0 F . It should be noted that for the 
orbit height under consideration, nadir angles greater than 70 deg. have not been considered as the 
line of sight does not contact the earth’s surface at larger angles. 
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Figure (1.10) Dependence of fore/aft optical axes angular separation on nadir angle. 



Nadir angle, 0 n 


Figure ( 1.11) Expansion of Figure ( 1 . 10) to show region of interest. 


Lula/ Performance Analysis 


37 






Figure (1.12) Dependence of the ground track length between the fore and aft line of sight 

ground points on nadir angle. 



Figure (1.13) Expansion of Figure ( 1 . 1 2) to show the region of interest. 
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Forward look time, (s) 



Figure (1.14) Dependence of forward look time on nadir angle. 



Figure (1.15) Expansion of Figure (1.14) to show the region of interest. 
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4.0 Appendix C - Lidar Performance Calculation 

4.1 Description 

This appendix lists the code used to generate the lidar model. Unlike conventional programming, 
spreadsheet programs do not consist of a linear list of commands but a series of expressions at 
various locations throughout the spreadsheet. The built-in code provided by the spreadsheet ven- 
dor enables a user iterface to be created with relative ease. The nature of this program which cen- 
ters on looking up values in assorted tables and then performing limited numerical calculations on 
them is ideally suited to a spreadsheet. Results can then easily be presented to the user through the 
use of the vendor supplied graph plotting routines. 

This program consists of three files :- 

• Lidar. WB 1 - the main spreadsheet which obtains user inputs, performs calculations and pre- 
sents results. 

• Atmos. WB 1 - a spreadsheet containing the atmospheric extinction as a function of altitude for 
a range of wavelengths. 

• Lidar.BAR - a toolbar which is used to provide access to the input routines. 

The program is written for the Borland Quattro-Pro for Windows spreadsheet. Each section of the 
code is contained on a seperate notebook page of the spreadsheet. The code on each page is listed 
below. A ‘user’s view’ of each page was provided in Section 1.9. 

4.2 Using the program 

To run the program, the user should place the three files listed above in a common directory. Quat- 
tro Pro for Windows should then be started and thefile Lidar. WB1 loaded. When loading this pro- 
gram will ask if the user wishes to open supporting hotlinks (this refers to the Atmos.WBl file). 
The answer to this question should be ‘OK’. The Lidar.BAR speedbar is loaded as a second 
speedbar by following the instructions in the Quattro Pro documentation. It is intended to auto- 
mate the start-up procedure at a later time. Once the program is loaded the user can select input 
parameters for the various lidar components by ‘clicking’ on the required component button on 
the speedbar and then changing the values in the dialog boxes. All plots update automatically 
(unless the user has explicitly turned off automatic update) except for the performance plot which 
is only updated when the ‘Update Graphs’ button on that page is selected. 

At the time of this report additional velocity estimation parameter regimes were still being pro- 
vided by Dr. Frehlich therefore a routine to automatically select the optimum parameterisation for 
a pair of ft and M values had not been implimented. This means that the user must manually tell 
the spreadsheet which parameterisation to use. This is done by typing the row number of the clos- 
est ftM pair on the velocity estimator look-up table into the cell labeled “MLE row no.” on the 
Results page of the spreadsheet. 

The code is still under development with current work concentrating on getting the code opera- 
tional rather than ‘bullet-proof’ and so none of the help buttons currently provide messages and 
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the user is not protected from modifying the spreadsheet. A backup copy of the files should there- 
fore be kept in another location. 

4.3 Variable list 

The following table lists the variables used and their locations. Spreadsheet cell references are of 
the form [SSNAME]PGNAME:CELL where SSNAME is the name of the relevant spreadsheet 
(defaults to active current spreadsheet if not supplied), PGNAME is the relevant page of the 
spreadsheet and CELL is the cell address on the page. 

The values in this table were exported from the spreadsheet page Block_Names. Due to differ- 
ences in the number of columns used in the spreadsheet and here the list is not in strict alphabeti- 
cal order, where alphabetical order is not followed, a double horizontal line occurs in the table 
below. 


Variable Name 

Cell Location 

Variable Name 

Cell Location 

ALIGNDAT 

Alignment:A36..B173 

EPSILON 

ML_Estimator:I2 

ALIGN_ANGLE 

Instrument:C22 

ETA_ ALIGN 

AIignment:B37 

ALIGN_LOSS 

Alignment:E44 

ETA_EXTTNCT 

Extinction:C242 

ALPHA 

ML_Estimator:E2 

ETA.HETERODYNE 

Instrument: G8 

ALT 

Results:C2 

ETA_HQE 

Instrument: G6 

ALTITUDE 

Extinction:A37 

ETA_MISALIGN 

Instruments 10 

ALTITUDE 1 

CoherenceL:A37 

ETA_POL ARISE 

Instrumental 

ATTEN 

Results:Cl3 

ETA.SYSTEM 

Instrument:G9 

B 

ResuIts:G5 

EXTINCTDAT 

Extinction: A36..G237 

BO 

ML_Estimator:D2 

FALSESPEC_MSG 

messages:D3..H3 

BETA 

Results:C3 

FREQ 

Instruments 14 

BETA_CALC 

Macros:B27 

G 

Constants :B7 

BETA_COUNT 

Macros:B9 

GO 

MLJEstimator: H2 

BETA_RANGE 

Perform:A36 ..A59 

GAMMA 

ML_Estimator:F2 

c 

Constants:B3 

GEOMETRY 

Instrument:GS 

CHI 

ML_Estimator:G2 

GEOMETRY.TYPES 

Lists :C2..G4 

CN2 

CoherenceL:E37 

H 

Constants:B2 

COHERENCE 

CoherenceL;C241 

HCOHERENCE 

Constants :B 8 
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D 

Instrument:C23 

LAMBDA 

Instrument:C8 

DELTA 

ML_Estimator:J2 

LAMBDA.LIST 

Lists: A1..A8 

DUMMY 1 

AIignment:E43 

LASER_DATA 

Instrument:C8..C12 

DUMMY2 

Conslants:Bl 1 

LATITUDE 

Scan_Pattcrn:038 

DUTY 

Instrument:Cl 1 

LOSRANGE 

Results:G2 

M 

Results:G6 

PHI 

Results :G7 

MARGIN 

Instrument:G24 

PLOTMINS 

Instrument:G18 

ME 

Constants:B5 

POLN 

Instrument:C21 

min_beta_point 

Macros:B10 

PRF 

Instrument:C12 

MIN_PROB 

Macros: B1 1 

QE 

Instrument:G6 

MLS.CLEAR 

Perform:B I..I202 

R 

Resuits:C10 

MSG1 

messages:B2..B7 

RANGE 

Extinction :D37 

MSG2 

messages:B9..B28 

RANGE2 

CoherenceL:C37 

MSG3 

messages:B30..B50 

RE 

Constants:B4 

MTS 

Rcsults:G2 

RECEIVER.DATA 

Instrument :G4 . .G6 

MU 

ML_Estimator:K2 

RECEIVER_GEOM 

Lists :C1..G4 

MULTIPLIER 

CoherenccL:F37 

RECEIVER_TYPE 

Lists:Il..I2 

NAD ALT 

Rcsults:C9 

RHOO 

Results:CU 

NADALT1 

Extinction:C37 

RO 

Instrument:C20 

NADALT2 

CoherenceL:B37 

RPM 

Instrument :G 19 

NADIR 

Instrument^ 18 

RSWATH 

Results:C26 

NAROW_BANDWIDTH 

Results:GI3 

SAVE_BETA 

Macros:B8 

NRG 

Instrumcnt:C9 

SCANNER.DATA 

InstrumentiG 1 7 . .G 1 9 

NS HOTS 

Results:Gl 1 

SCAN.DATA 

Instrument:G17..G19 

ORBH 

Instrument:C4 

SCAN.FLAG 

Scan_Pattem:J39 

ORBINC 

Instnjment:C5 

SCAN.TYPE 

Instrument:G17 

ORBIT_DATA 

Instrument:C4 ,C5 

SCAN_TYPE_LIST 

Lists J1.J3 

PARAM 

Alignment:A37 

SIGMAM 

Results:G24 

PATH 

Extinction:E37 

SIGMAS 

Results:C7 














































































































PATH2 

CoherenceL;D37 

SIGMA VLOS 

Results:G25 

PBAD 

Results:G22 

SIGMA V_H 

Results:C5 

PERFORMANCE 

Perform:A36..E59 

SIGMA V_I 

Results:G24 

PERFORM_CALC 

Macros :B 13 

SIGMAV.V 

Results:C6 

PGOOD 

ResuIts:G23 

SLANTG 

Results:C23 

SNR 

Results:G20 

TSTEP 

Scan_Pattem:B 1 

SNR_NARROW 

Resuits:G21 

TT 

Instrument:G7 

SNR.WIDE 

Results:G20 

TTAR 

Instrument:C19 

SPOTG 

Results:C24 

T_S PEC_WIDTH 

Constants :B9 

SYSTEM 

Instrument :G9 

VI 

Results:C21 

TARGET.ALTTrUDE 

Results:C2 

VEARTH 

Results:C22 

TARGET.DATA 

Results:C2..C8 

VERT.RANGE 

Results:C8 

TAU 

Instrument:C10 

VG 

Results :C21 

TE 

Constants:B6 

VMAXH 

Results:C4 

TEST 

Extinction:L27 

VMAXLOS 

Results :C 14 

TFLIP 

Results:C27 

VSAT 

Results:C20 

TO 

Instrument:Cl9 

VSEARCH 

Results :C4 

TORBIT 

Results:C25 

w 

Results :G8 

TRAR 

Instrument:C20 

WIDE.BANDWIDTH 

Results:G5 

TS 

Results:G4 

WPL37DAT 

CoherenceL:A36..F237 

TSCOPE_D 

Instrument:C17 

\P 

Perform: 140 


4.4 Instrument 


The following code lists the contents of cells on the instrument page. Text cells contain the initial 
characters, \ A and “ which indicate left, center and right alignment respectively. Cell contents 
beginning wtih a + or @ are calculation cells. All of the cells are visible to the user. 


Instrument : Cl 
Instrument : B3 
Instrument : F3 
Instrument : B4 
Instrument : C4 
Instrument : D4 
Instrument : F4 
Instrument : G4 


' Instrument Parameters 

"ORBIT 

"RECEIVER 

"Orbit height 

350 

' km 

"Type 

' Complex 
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Instrument : B5 : "Inclination angle 

Instrument :C5 : 98 
Ins trument : D5 : ' deg 

Instrument : F5 : "Geometry 
Instrument :G5 : 'Wang 

Instrument : F6 : "Detector quantum efficiency 
Instrument : G6 : 0.6 
Instrument : B7 : ''LASER 

Instrument : F7 : "Transmit beam intensity fraction 

Instrument : G7 : (3VLOOKUP ( $GEOMETRY , $GEOMETRY_TYPES , 1 ) 

Instrument : B8 : "Wavelength 
Instrument : C 8 : 2.065479 
Instrument : D8 : 'urn 

Instrument : F8 : "Heterodyne efficiency 

Instrument : G8 : @VLOOKUP ( GEOMETRY , GEOMETRY_TYPES , 2 ) 

Instrument : B9 : "Pulse energy 

Instrument : C9 : 0.025 
Instrument : D9 : 'J 

Instrument : F9 : "System efficiency 
Instrument : G9 : +ETA_HETERODYNE*TT 
Instrument : BIO : "Pulse length 

Instrument : CIO : 0.5 
Instrument : D10 : 'us 

Instrument : F10 : "Misalignment efficiency 
Instrument :G10 : +ALIGN_LOSS 
Instrument : Bll : "Duty cycle 
Instrument :C11 : 0.05 
Ins trument : B1 2 : "P.R.F. 

Instrument : C12 : 10 

Instrument : D12 : 'Hz 

Instrument : B13 : "Spectral width 

Instrument :C13 : +T_SPEC_WIDTH 

Ins trument : D13 : 'MHz 

Instrument : B14 : "Frequency 

Ins trument : C14 : +C/LAMBDA 

Instrument : D14 : 'MHz 

Instrument : B16 : A OPTICS 

Instrument : F16 : “'SCANNING 

Instrument : B17 : "Telescope diameter 

Instrument : C17 : 0.25 

Instrument : D17 : 'm 

Ins trument : F17 : "Scan type 

Instrument : G17 : 'Wedge 

Ins trument : B1 8 : "Nadir angle 

Instrument : C18 : 30 

Ins trument : D1 8 : 'deg 

Ins trument : FI 8 : "Plot duration 

Instrument : G18 : 1 

Instrument : H18 : 'mins 

Ins trument : B19 : "Transmit optics 

Ins trument : C19 : 0.9 

Instrument : F19 : "Telescope rotation rate 

Instrument :G19 : 12 

Instrument : H19 : ' rpm 

Instrument : B20 : "Receive optics 

Ins trument : C2 0 : 0.9 

Instrument : F20 : "{ Conical /wedge scan only) 

Instrument : B21 : "Polarisation efficiency 

Ins trument : C2 1 : 1 

Instrument : B22 : "Receive/ lo misalignment angle 

Inst rument : C22 : 7.2819023506228 

Instrument : D22 : 'prad 

Instrument : B23 : "Effective diameter 

Instrument :C23 : @ I F ( SCAN_TYPE= * Wedge - , +TSCOPE_D*9COS (©RADIANS ( $NADIR) ) , TSCO 
PE_D) 
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Instrument : D23 
Instrument : F23 
Instrument : B24 
Instrument : F24 
Instrument : G24 
Instrument : B2 5 
Instrument : C25 
Instrument : D2 5 
Instrument : B26 
Instrument : C2 6 
Instrument : D2 6 


' m 

■"SYSTEM 

"(smaller for wedge scan) 

"System margin 

0.5 

"Rayleigh range 
+ D A 2 / ( @ PI * LAMBDA* 1 0 A -3 ) 

' km 

"Beam spot size 
+D/@PI 
' m 


4.5 Results 


The following code lists the cell contents of the results page. All of the cells are visible to the user. 

Results : B1 : A TARGET 
Results: FI: A SIGNAL PROCESSING 
Results :B2: "Aerosol altitude 

Results :C2 : 300 
Results : D2 : 1 m 

Results :F2: "Line of sight range resolution 

Results :G2 : +VERT_RANGE/ (2COS { @ RADIANS (NADALT) ) 

Results :H2*: 'm 

Results :B3: "backscatter (lambda) 

Results:C3: 4.9E-06 
Results:D3: '/ (m-sr) 

Results :F3: "Observation time 

Results :G3 : + L0SRANGE* 2 / (C*10 A -6) 

Results : H3 : 'us 

Results :B4: "Maximum horizontal velocity 

Results : C4 : 30 
Results : D4 : 1 m/ s 

Results :F4: "Time between sample points 

Results :G4 : + LAMBDA/ VMAXLOS *@IF ( Instrument : $.G$ 4 = "Complex* ,4,8) 

Results : H4 : 'us 

Results :B5: "Horizontal wind velocity uncertainty 

Results : C5 : 0 
Results : D5 : ' m/ s 

Results :F5: "Digitisation frequency 

Results :G5: 1/TS 

Resu 1 t s : H5 : 'MHz 

Results :B6: "Vertical wind velocity uncertainty 

Results:C6: 0 
Results : D6 : ‘ m/ s 

Results :F6: "No. sample points / observatn . 

Results:G6: +G3/TS 

Results :B7: "Wind variance between shots 

Results:C7: 0 
Resul ts : D7 : * m/ s 

Results : F7 : " Phi 

Results :G7 : 10 A (SNR/ 10) *M/@IF < Instrument : $G$4= "Complex* ,1,2) 

Results :B8: "Vertical range resolution 

Results : C8 : 1000 

Results : D8 : * m 

Results : F8 : "Signal width 

Results :G8 : @SQRT ( (0.1874/ (TAU) ) A 2 + { ( SIGMAV_H/0SIN ( SRADIANS (NADALT) ) ) *2/ 
LAMBDA) A 2 + ( { S IGMAV_V/ (iCOS ( @ RAD IANS ( NADALT ) ) ) * 2 / LAMBDA ) A 2) 

Results : H 8 : ' MHz 

Results :B9: "Target nadir angle 

Resul t s : C9 : ^DEGREES ( £AS IN { (ORBH* 10~3 + RE) * @S IN ( ^RADIANS (NADIR) ) / (RE+ALT) ) ) 

Resul t s : D9 : ' deg 

Results :F9: "Omega 
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Results :G9: +W*G3 
Results:H9: 'm/s 

Results: BIO: "Line of sight range to this altitude 

Results : CIO : 0SQRT ( < $ORBH*10"3+$RE) "2+{$RE+ALT) "2- ( 2 *$ORBH*10"3+2 *$RE) * 

{ $RE+ALT) *@COS ( @ RADIANS (NADALT-$NADIR) ) ) /1000 
Results: DIO: 'km 

Results : F10 : "Sigmav/w 

Results :G10 : @ ARRAY ( (+@@ ( "ML_Es timator : G* &<£STRING (G14 , 0) ) * ( 1+ ( PHI/ 

( "ML_Es timator : H"&@STRING (G14 , 0) ) ) "+@@ ( "ML_Es timator : I *&@STRING (G14 , 0) ) ) " ( - 
@@ ( "ML_Es timator : J*&@STRING (G14 , 0 ) ) ) + + { "ML_Estimator : K"&(iSTRING (G14 f 0 } ) ) } 

Results :B11: "Coherence length 

Results :C11: +COHERENCE 
Results :D11: 'm 

Results :F11: "No. of shots/ wind estimate 

ResultsrGll: 1 

Results :F12: "Bandwidth (wide band) 

Results :G12 : 4 *VMAX LOS/ LAMBDA 
Results : H12 : 'MHz 

Results :B13: "One way Intensity Transmission 

Results : C13 : +ETA_EXTINCT 

Results :F13: "Bandwidth (narrow band) 

Results : G13 : +W 
Results : H13 : ' MHz 

Results :B14: "Maximum line of sight velocity 

Results :C14 : +VSEARCH* @SIN ( @ RADIANS ( NADALT) ) 

Results:D14: 'm/s 

Results :F14: "MLE row no. 

Results :G14: 6 

Results :B19: "OTHER PARAMETERS 
Results : F19 : "RESULTS 
Results :B20: "Satellite velocity 
Results :C2 0 : @SQRT (G*ME/ (ORBH* 10 " 3 +RE) ) 

Results:D20: 'm/s 

Results :F20: "Wideband SNR 

Results :G20 : 10*@LOG ( ( NRG*TT*TO*ATTEN"2 *BETA*@PI *D"2 *RO* SYSTEM *ETA_MI SAL IGN* 
POLN*QE*C*$MARGIN) / (8*H* (FREQ*10"6) *$G$12*10"6* ( R* 1000 ) "2+ ( 1+ ( D ( 2 *RHO0 ) ) "2) ) ) 
Resul t s : H2 0 : ' dB 

Results :B21: "ground track velocity 

Results :C21 : +VSAT*RE/ ( ORBH * 1 0 " 3 +RE ) 

Results:D21: 'm/s 

Results :F21: "Narrowband SNR 

Results :G21 : 10*@LOG{ 10" (SNR/ 10) *2*VMAXLOS/ ( LAMBDA* W) ) 

Results : H2 1 : ' dB 

Results :B22: "Earth rotation velocity at equator 

Results :C22 : 2*(iPI*RE/TE 
Results : D2 2 : 'm/s 

Resul t s : F2 2 : "P(bad) 

Results :G2 2 : @ ARRAY ( ( 1+ ( PHI / ( "ML_Es timator : D"&@STRING (G14 , 0 ) ) ) "ML_Es- 
t imator : E"&@ STRING (G14 , 0 ) ) ) "- + @@ ( "ML_Es timator : F"$c@STRING (G14 , 0 ) ) ) 

Results :B23: "Slant range to ground 

Results :C2 3 : @SQRT( (ORBH* 10"3+RE) "2 + RE"2- ( 2 *ORBH* 10"3 + 2 *RE) *RE*3COS ( GRADIAN- 
S (NADALT -NADIR) ) ) 

Results : D23 : 'm 

Results : F2 3 : "P (good) 

Results :G23 : 1-PBAD 

Results :B24: "Spot diameter on ground 

Results :C2 4 : + $ D * @ SQRT { 1 + ( LAMBDA * 1 0 " - 6 * SLANTG / (0PI* ($D/2) "2) ) "2) 

Results : D2 4 : ’ m 

Results :F24: "sigmavlos - instrument 

Results : G 2 4 : +G1 0 *W* LAMBDA/ 2 

Results : H2 4 : ' m/ s 

Results :B25: "Time for one orbit 

Results : C 2 5 : 2* (2 PI* (ORBH* 10 " 3 *RE)./ VSAT 

Resul ts : D2 5 : ’ s 

Results : F25 : "sigmavlos 
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Results :G 25 
Results : H25 
Results : B2 6 
Results : C2 6 
Results : D26 
Results : F26 
Results : G26 
Results : B27 
Results : C27 
Results : D27 
Results : F27 
Results : G27 


@SQRT ( { SIGMAV_I /s 2 + SIGMAS / '2 ) / (NSHOTS* < 1-PBAD) ) ) 

'm/s 

"Swath radius (conical /wedge scan) 

2 *RE*@RADIANS (NADALT- NADIR) / 2 
'm 

"F.O.M. 

( (ATTEND ) / (R A 2) ) * <*S IN ( (iRADIANS (NADALT) ) *10^6 
"Optimum mirror flip time (line scan) 

2 *RE/V1*@ASIN ( SLANTG*@SIN ( @ RAD IANS (NADIR) ) / (@SQRT(2) *RE) ) 
's 

"PHI 
+ PHI 


4.6 Scan pattern 


The following code lists the contents of the Scan_Pattem page. The user sees two plots, one of the 
scan pattern and the other of the beam separation for the line scan pattern. 


The following lines appear above the plots. 


Scan_Pattern : A1 
Scan_Pattern : B1 
Scan_Pattern : Cl 
Scan_Pattern : D1 
Scan_Pattern : II 


' tstep 

1/ ( PLOTMINS * 60 / 600 ) 

' Hz 

'Shot Pattern at Equator 

'Fore/Aft Look Seperation Dependence on Latitude 


The following line is visible to the user below the plots. 

Scan_Pat tern : 13 3 : ' This plot is relevant to the half of 'BEST' shot pattern. 


The remaining lines are out of sight from the user and are used to calculate the scan pattern. The 
model assumes that the satellite starts from 0 deg. latitude and O.deg longitude. Columns B and C 
calculate the spacecraft latitude and longitude, column D calculates the telescope pointing angle 
with respect to the spacecraft velocity vector and, columns E and F calculate the position of the 
lidar beam on the ground (in km.) from the 0 deg. latitude, 0 deg. longitude position. Columns G 
and H calculate the beam separation for the line-scan pattern as a function of latitude. A total of 
600 shot positions are calculated. 


Scan_Pa ttern : B3 6 : 
Scan_Pat tern : D3 6 
Scan_Pat tern : B37 
Scan_Pattern : C3 7 


satellite track 
' angle 
' long 
' lat 


Scan_Pat tern : E37 
Sc an_Pat tern : F37 
Scan_Pa ttern : G37 


' long 
' lat 
' lat 


Scan_Pattern : H37 


1 separation 


Scan_Pat tern : A3 8 : 0 

Scan.Pa ttern: B3 8 : + $A3 8 /TSTEP* $ VI * @COS ( @ RADIANS ( 180-$ORBINC) ) /1000 + $A3 8/ 

$ TSTEP* SVEARTH * @COS ( @ RADIANS (SLATITUDE) ) /1000 
Sc an_Pa ttern : C3 8 : ♦ $A3 8 /TSTEP *$ VI *@SIN ( ©RADIANS ( 180 -SORBINC ) ) /1000 

Sea n^Pa ttern : D3 8 : @IF ( @MOD ( $A3 8 / STSTEP , STORBIT) > STORBIT* $ DUTY , 0 , GIF ( $SCAN_- 
FLAG = 1 , $RPM/ 60*2*@PI*$A38/ STSTEP+0 , 0 IF ( @MOD ( @INT ( $A3 8 / $TSTEP/ $TFLIP) , 2) >0 , £ RA- 


DIANS (13 5) , 0 RADIANS (45) } ) ) 

Scan_Pat tern : E38 : + $B3 8 ♦ $RSWATH*@COS ( @PI -(2RADIANS ( SORBINC) +SD38) /1000 

Scan_ Pa ttern: F3 8 : + $C3 8 + SRSWATH* @SIN ( (iPI -@RADIANS ( SORBINC ) +$D3 8 ) /1000 


Scan_Pat tern : G3 8 : 0 

Scan_Pat tern : H3 8 : + SVEARTH * $TFLIP*@COS ( @ RADIANS ( $G3 8 ) ) /1000 
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Scan_Pattern : N3 8 : 'latitude 

Scan_Pattern:038 : 0 
Scan_Pattern: P38 : 'deg 

Scan_Pattern : A3 9 : 1 

Scan_Pattern : B3 9 : +$A39 /TSTEP*$V1 *@COS ( ©RADIANS ( 180-$ORBINC) ) /1000+SA39/ 
$TSTEP*$VEARTH*0COS (0RADIANS ( $LATITUDE) ) /1000 
Scan_Pat tern : C3 9 : +$A39/TSTEP*$V1*@SIN(@RADIANS(180-$ORBINC) ) /1000 
Scan_Pattern : D3 9 : ©IF ( ©MOD ( $A39 /STSTEP , STORBIT) >$TORBIT*$DUTY,0,@IF($SCAN_- 
FLAG=1, $RPM/60*2*@PI*$A39/$TSTEP+0, @IF(@MOD(@INT($A39/$TSTEP/$TFLIP) , 2) >0, ©RA- 
DIANS (135), ©RADIANS ( 45 ) ) ) ) 

Scan_Pattern : E3 9 : +$B3 9 +SRSWATH* ©COS ( ©PI -©RADIANS ( SORBINC) +$D3 9) /1000 
Scan_Pattern:F39 : +$C39+$RSWATH*@SIN ( ©PI-0RADIANS ( SORBINC) +SD39) /1000 
Scan_Pattern:G39 : 5 

Scan_Pattern:H39 : +SVEARTH* $TFLIP*@COS ( ©RADIANS ( $G3 9 ) ) /1000 
Scan_Pattern: J39 : ©IF ( SCAN_TYPE=" Dual- look" ,0,1) 

Scan_Pattern: A40 : 2 

Scan_Pattern : B40 : + $A4 0 / TSTEP *$ VI* ©COS ( ©RADIANS ( 180-$ORBINC) ) /1000+SA40/ 
$TSTEP*$VEARTH*8COS (©RADIANS (SLATITUDE) ) /1000 
Scan_Pattern : C40 : +$A40/TSTEP*$V1*@SIN(@RADIANS (180-SORBINC) ) /1000 
Scan_Pattern : D40 : ©IF ( ©MOD ( SA4 0 / $TSTEP , $TORBIT) >$TORBIT*$DUTY , 0 , ©IF ( $SCAN_- 
FLAG=1 , $RPM/ 60*2 *@PI* $A40 / STSTEP+0 , ©IF ( ©MOD ( ©INT ( $A4 0/ STSTEP / STFLIP) , 2) >0, ©RA- 
DIANS (13 5) , ©RADIANS (45) ) ) ) 

Scan_Pa t tern :E40: + $ B4 0 + $ RSWATH * ©COS ( © PI - ©RADIANS ( $ORBINC )+$D40)/1000 
Scan_Pattern: F40 : +$C40+$RSWATH*@SIN(@PI-©RADIANS ($ORBINC) +$D40) /1000 
Scan_Pattern : G40 : 10 

Scan_Pattern : H40 : +SVEARTH* STFLIP* ©COS ( ©RADIANS ( $G40 ) ) /1000 


The lines A40 - H40 are repeated, with column A being incremented by 1 each row and column G 
being incremented by 5 each row down to the following rows:- 

Scan_Pattern : A55 : 17 

Scan_Pattern:B55 : + $A55 /TSTEP* $Vl*@COS (©RADIANS ( 180-$ORBINC) ) /1000+SA55/ 
STSTEP *$VEARTH* ©COS ( ©RADIANS ( $LATITUDE) ) /1000 
Scan_Pattern:C55 : +$A55/TSTEP*$V1*@SIN(@RADIANS(180-$ORBINC) ) /1000 
Scan_Pattern : D55 : ©IF ( ©MOD ( $A55 / STSTEP , $TORBIT) > STORBIT* $DUTY, 0 , ©IF ( $SCAN_F 
LAG=1 , $RPM/ 60*2*@PI* $A55 / STSTEP+0 , ©IF ( ©MOD (0 INT ( $A55/$TSTEP/$TFLIP) , 2) >0, ©RA- 
DIANS (13 5) , ©RADIANS (45) ) ) ) 

Scan_Pattern : E55 : +$B55+$RSWATH*@COS ( ©PI -©RADIANS ( SORBINC) +$D55) /1000 
Scan_Pattern: F55: + SC55 + $RSWATH*@SIN ( ©PI -©RADIANS ( SORBINC) +SD55) /1000 
Scan_Pattern : G55 : 85 

Scan_Pattern:H55 : + SVEARTH* STFLIP* ©COS (©RADIANS ($G55) ) /1000 
Scan_Pattern: A56 : 18 

Scan_Pattern: B56 : + $A56 /TSTEP* $V1 * ©COS ( ©RADIANS ( 180-SORBINC) ) / 1000-*-$A56/ 

$TSTEP* $VEARTH * ©COS ( ©RADIANS ( $ LATITUDE ) ) /1000 
Scan_Pattern : C56 : +$A56 /TSTEP*$V1 *©SIN ( ©RADIANS ( 180 -SORBINC) ) /1000 
Scan_Pattern : D56 : ©IF ( @MOD( $A56/ $TSTEP , STORBIT) >$TORBIT*$DUTY, 0 , ©IF ( $SCAN_- 
FLAG=1 , $RPM/60*2*@PI*$A56/$TSTEP+0, ©IF ( ©MOD ( ©INT ( $A56/$TSTEP/ $TFLIP) , 2) >0, ©RA- 
DIANS (13 5) , ©RADIANS (45) ) ) ) 

Scan_Pattern : E56 : +$B56+$RSWATH* ©COS ( ©PI -©RADIANS ( SORBINC) +$D56) /1000 
Scan_Pattern:F56 : +$C56 +SRSWATH* ©SIN ( ©PI -©RADIANS ( SORBINC) +SD56) /1000 

Scan_Pattern : G56 : 90 

Scan_Pattern:H56 : +SVEARTH* STFLI P*@COS ( ©RADIANS ( SG56 ) ) /1000 


This completes the line-scan beam separation calculation and from this point only the shot pattern 
calculation is repeated until:- 

5can_Pattern : A637 : 599 

3can_Pattern : B637 : + $A637 /TSTEP* $V1 *@COS ( ©RADIANS ( 180 -SORBINC) ) /1000 + SA637/ 

STSTEP 'SVEARTH* ©COS (©RADIANS ( SLATITUDE) ) /1000 
3can_Pat tern : C6 3 7 : + $A637 / TSTEP* $V1 *©SIN ( ©RADIANS ( 180-SORBINC ) ) /1000 
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Scan_Pattern : D637 : GIF ( GMOD ( $A637 / $TSTEP , $TORBIT) >$TORBIT*$DUTY , 0 , GIF < $SCAN_- 
FLAG=1 , $RPM/ 60 *2 *GPI * $A637/ $TSTEP + 0 , GIF { GMOD ( GINT ( $A63 7 / $TSTEP/ $TFLIP ) , 2) >0, G- 
RADIANS (13 5) f GRADIANS ( 45 ) ) ) ) 

Scan_Pattern : E637 : +$B637+$RSWATH*GCOS ( GPI-GRADIANS ( $ORBINC) +$D637) /1000 
Scan_Pattern:F637 : +$C637+$RSWATH*GSIN ( GPI-GRADIANS ( $ORBINC ) +$D637 ) /1000 
Scan_Pattern : A638 : 600 

Scan_Pattern:B638 : +$A638/TSTEP*$V1*GC0S (GRADIANS (180-$ORBINC) ) /1000+$A638/ 
$TSTEP* $VEARTH*GCOS ( GRADIANS { $LATITUDE) ) /1000 
Scan_Pattern: C63 8 : +$A638/TSTEP*$V1*GSIN ( ©RADIANS ( 180-$ORBINC) ) /1000 
Scan_Pat tern : D63 8 : GIF ( GMOD ( $A63 8 / $TSTEP , $TORBIT) >$TORBIT* $DUTY, 0 f GIF ( $SCAN_- 
FLAG=1 , $RPM/ 60 * 2 *GPI * $A63 8 / $TSTEP+0 , GIF (GMOD(GINT($A638/$TSTEP/$TFLIP) , 2) >0, G- 
RADIANS (135) , GRADIANS (45) ) ) ) 

Scan_Pattern:E638 : + $B63 8 + $RSWATH*GCOS ( GPI - GRADIANS ( $ORBINC ) +$D638) /1000 
Scan_Pattern : F63 8 : + $C63 8+ $RSWATH*GS IN (GPI-GRADIANS ($ORBINC) +$D638) /10 00 


4.7 Extinction 


The following is a list of the code on the Extinction page. The user sees a plot of the one-way and 
two-way transmission of the atmosphere over the altitude range 0-20 km. The following lines are 
visible to the user beneath the plot. 

Extinction: D33 : 'Optical transmission from satellite to target, laser wave- 

length is 

Extinction : K3 3 : + LAMBDA 

Ext inction : L33 : 'um 


The following lines are not visible to the user and are used to calculate the the one-way and two- 
way atmospheric extinction. Column A contains the altitude above the earths surface (incre- 
mented in 0. 1 km steps), column C containd the nadir angle at that altitude (accounting for the 
earths curvature), column D contains the range from the satellite to the altitude and column E the 
path length from one altitude increment to the next. Column F contains the atmospheric extinction 
at the current altitude, column G the product of the path length and extinction, column H the sum- 
mation of the path length/extinction products from the spacecraft to the current altitude. Column I 
contains the one-way extinction to the altitude and column J the two-way extinction. 


Extinction : A3 6 : 
Extinction : B3 6 : 
Extinction : C3 6 : 
Ext inction : D3 6 : 
Extinction : E3 6 : 
Extinction : F3 6 : 
Extinction : G36 : 
Extinction : 136 : 
Extinction : J36 : 
Ext inction : A37 : 
Extinction : B37 : 
Extinction:C37 : 
TUDE *1000) ) ) 
Extinction : D37 : 
( 2 * SORBH *10 A 3>2* 
Extinction : E3 7 : 
Extinction : F37 *. 
Extinction : G3 7 : 
Extinction : H3 7 : 
Extinction : 137 : 
Ext inction : J3 7 : 
Extinction : A3 8 : 
Extinction : B3 8 : 


' Altitude 
' Index 
’ Nadaltl 
' range 
' path 
+ LAMBDA 
"lamda*pathl 
'One Way 
'Two Way 
0 
1 

GDEGREES { GASIN ( ( ORBH* 10 A 3 +RE) *GSIN ( GRADIANS (NADIR) ) / (RE+ALTI- 

GSQRT ( ( SORBH* 1CP3+SRE) ~2> ( $RE+ALTITUDE* 1000 ) ~2- 

$RE ) * ( $ RE + ALT I TUDE* 1000 ) *GCOS ( GRADIANS (NADALT1 -$NADIR) ) ) /1000 

0 

(iHLOOKUP ( $ LAMBDA , ( \QPW\LIDAR\ ATMOS] $MLS_CLEAR, $B37) 

0 

GSUM { $G3 7 . . $G$2 3 7 ) 

GEXP { -H37 ) 

>137^2 
0 . 1 
2 
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Extinction : C3 8 : ©DEGREES { ©AS IN ( { ORBH* 1 0 A 3 +RE ) * ©S IN ( ©RADIANS (NADIR) ) / 
(RE+$A38*1000 ) ) ) 


Extinction :D3 8 : ©SQRT ( ( $ORBH* 10^ 3 + $RE) "2 + ( $RE+$A38* 1000 ) "2- ( 2 *$ORBH* 10^3 + 2* 
$RE) * ( $RE+$A3 8*1000) * ©COS ( ©RADIANS ( $C3 8- $NADIR) ) ) /1000 


Extinction:E38 
Extinction : F3 8 
Extinction : G3 8 
Extinction : H3 8 
Extinction : 1 3 8 
Extinction : J3 8 


+RANGE-$D38 

©HLOOKUP ( $ LAMBDA, [ \QPW\LIDAR\ATMOS ] $MLS_CLEAR, $B38) 

+ $E3 8 * $F3 8 

©SUM ( $G3 8 . . $G$237 ) 

©EXP ( -H3 8 ) 

+138^2 


A38 - J38 are repeated (withsuitable incrementing) until 


Extinction : A23 6 : 
Extinction : B23 6 : 
Extinction : C23 6 : 
( RE+$A2 3 6*1000) ) ) 
Extinction :D2 3 6 : 
2 * $RE ) * ( $RE + $A2 3 6 
Extinction: E236 : 
Extinction : F236 : 
Extinction : G23 6 : 
Extinction : H2 3 6 : 
Extinction : 1236 : 
Extinction'': J23 6 : 
Extinction : A2 3 7 : 
Extinction : B237 : 
Extinction : C237 : 
( RE + SA237 * 1000 ) ) ) 
Extinction : D2 37 : 
2 * $RE ) * ( $RE+ $A2 3 7 
Extinction : E2 3 7 : 
Extinction : F237 : 
Extinct ion : G237 : 
Ext inct ion : H2 3 7 : 
Extinction : 12 3 7 : 
Ext inct ion : J2 3 7 : 


19 . 9 
200 

©DEGREES (©ASIN{ ( ORBH* 10 A 3 + RE) * ©SIN { ©RADIANS (NADIR) ) / 

©SQRT ( ( $ORBH* 10^3 + $RE ) ^2+ ( $RE+$A236*1000 ) A 2 - ( 2 * $ORBH* 10^ 3 + 
1000} *©COS ( ©RADIANS ( $C236-$NADIR) ) ) /1000 
+ $D23 5-$D23 6 

©HLOOKUP ( $LAMBDA, [ \QPW\L IDAR\ ATMOS] $MLS_CLEAR, $B236) 

+$E236 *$F236 

©SUM ( $G2 3 6 . . $G$23 7 ) 

©EXP ( -H2 3 6 ) 

+1236^2 

20 

201 

©DEGREES (@ASIN( ( ORBH* 10 A 3 + RE) * @S IN ( ©RADIANS (NADIR) ) / 

©SQRT ( ( $ORBH* 10 ^3 + $RE ) A 2 + ( $RE+$A237 * 1000 ) A 2- ( 2 *$ORBH*10 / '3 + 
1000) *©COS ( ©RADIANS ( $C237-$NADIR) ) ) /1000 
+ $D2 3 6 - $D237 

©HLOOKUP ($LAMBDA, [ \QPW\LIDAR\ ATMOS ] $MLS_CLEAR, $B237) 

+$E237 * $F237 

©SUM ( $G2 3 7 . . $G$2 3 7 ) 

©EXP ( -H237 ) 

+1237^2 


The following lines set up a search criteria and locate the extinction value for the altitude of inter- 
est (entered on the Results page). 

Extinct ion : A24 1 : ‘Search Criteria 

Ext inct ion : C2 4 1 : 'Result 

Ext inct ion : A2 42 : 'Altitude 

Extinction :C2 4 2 : ©EXP ( -©DSUM ( EXTINCTDAT , 6 , A242 . .A243) ) 

Extinction: A243 : +ALTITUDE> ( $ALT/ 1000) 

Extinction : C245 : 'This page searches the extinction database to find the alti- 

tudes that are greater 

Ext inct ion : C246 : 'than the target altitude. The values of extinction for the 

path length through each 

Ext inct ion : C247 : '0.1 km altitude band up to 20 km are summed. The value 

obtained is the value of 

Extinction :C248 : ' (alpha*path length) which is used to calculate the one way 

extinction through 

Ext inct ion : C249 : 'the atmosphere. 
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4.8 Coherence length 


The folowing code lists the contents of the cells on the page CoherenceL. The user sees a plot of 
the WPL-37 C n 2 profile and a plot of the transverse coherence length as a function of altitude 
from the earths surface. Column A contains the altitude, column B the nadir angle at that altitude, 
column C the range from the satellite to that altitude and column D the path length through an alti- 
tude increment. Column E calculates the refractive index structure function, C„ at the altitude 
and columns F and G calculate the transverse coherence length. 


CoherenceL : A3 6 
CoherenceL : B3 6 
CoherenceL : C3 6 
CoherenceL : D3 6 
CoherenceL : E3 6 


' Altitudel 
' Nadalt2 
’ range2 
' path2 
' Cn2 


CoherenceL : F3 6 
CoherenceL : A3 7 
CoherenceL : B3 7 
TUDE1*1000 ) ) ) 


'multiplier 

0 

©DEGREES ( ©ASIN ( ( ORBH* 1 0 " 3 +RE) *0SIN(0RADIANS (NADIR) ) / (RE+ALTI- 


CoherenceL : C37 : ©SQRT ( ($ORBH*10"3+$RE) "2 + ( SRE+ALTITUDEl* 1000 ) "2- 
(2*$ORBH*10"3+ 2*$RE) * < $RE+ALTITUDE1*1000) * ©COS (©RADIANS (NADALT2 -$NADIR) ) ) / 
1000 


CoherenceL : D37 : 0 

CoherenceL: E3 7 : ©IF ( ALTITUDE1>20 , 0 , ©IF ( ALTITUDE1<0 .01, 7 .36*10" (-14) , 0IF(ALTI- 
TUDE1<3 ,1.59 *10"- 12* (ALTITUDE1*1000) " (-4/3) , 3 . 66*10" (-17) ) ) ) 

CoherenceL: F3 7 : ©IF ( RANGE2>$R, 0 , PATH2 *10 "3 *CN2 * ( 1-RANGE2/$R) " (5/3) ) 

CoherenceL :G3 7 : ( $HCOHERENCE* (2*@PI/ ( $LAMBDA*10" ( -6 ) ) ) "2* 

©SUM ( F37 . . $F$237 ) )"(-3/5) 

CoherenceL : A3 8 : 0.1 

CoherenceL : B38 : ©DEGREES ( ©ASIN ( (ORBH* 10 "3 +RE) *0SIN ( 0RADIANS (NADIR) ) / 
(RE+$A38*1000) ) ) 

CoherenceL :C3 8 : 0SQRT( ( $ORBH*10"3 +$RE) "2+ ($RE+$A38*1000) "2- (2*$ORBH*10"3+ 
2*$RE) * ( $RE + $A3 8*1000) *0COS ( 0RADIANS ( $B3 8-$NADIR) ) ) /1000 
CoherenceL :D3 8 : +RANGE2 - $C3 8 

CoherenceL : E3 8 : 0IF ( $A3 8>20 , 0 , 0IF ( $A3 8<0 . 01 , 7 .36*10" ( -14) , 

0IF ( $A38<3 , 1 . 59*10"-12* ( $A3 8* 1000 ) " ( -4/ 3 ) , 3 . 66*10" (-17) ) ) ) 

CoherenceL :F3 8 : 0IF ( $C38>$R, 0 , $D38*10"3 *$E38* ( 1-$C38/$R) " (5/3 ) ) 

CoherenceL :G3 8 : ( SHCOHERENCE* (2*0PI/ ( $ LAMBDA* 10" ( -6) ) ) "2* 

@SUM(F38 . . $F$237 ) ) " ( -3/5) 

CoherenceL : A3 9 : 0.2 

CoherenceL : B39 : ©DEGREES ( 0ASIN ( (ORBH*10"3+RE) * §S IN (0RADIANS (NADIR) ) / 

(RE+SA39 * 1000 ) ) ) 

CoherenceL :C3 9 : ©SQRT ( ( $ORBH* 10 "3 + $RE ) "2 + ( $RE+$A39*1000 ) "2- ( 2*$ORBH*10"3 +2 * 
$RE ) * ( $RE+ SA39*1000) * ©COS ( ©RADIANS ( $B3 9- SNADIR) ) ) /1000 
CoherenceL : D3 9 : -t-$C38-$C39 

CoherenceL: E3 9 : ©IF ( $A3 9>20 , 0 , 0IF ( $A3 9<0 . 01 . 7 . 36*10" ( -14) , 

©IF ( $A39<3 , 1 . 59*10"-12* ( $A3 9*1000) " (—4/3 ) , 3 . 66 * 10" ( -17 ) ) ) ) 

CoherenceL: F3 9 : ©IF ( $C3 9>$R, 0 , $D39* 10 " 3 * $E39 * ( 1-$C39/$R) " (5/3) ) 

CoherenceL :G39 : (SHCOHERENCE* (2*@PI/ ( SLAMBDA* 10" ( -6) ) )"2* 

©SUM ( F39 . . $F$237 ) ) " (-3/5) 


A39 - G39 are repeated (with suitable increments) until:- 

CoherenceL : A236 : 19.9 

CoherenceL : B23 6 : ©DEGREES ( ©ASIN ( (ORBH*10"3+RE) *0SIN ( 0RADIANS (NADIR) ) / 

( RE + SA236 * 1000 ) ) ) 

CoherenceL :C2 3 6 : ©SQRT ( ( SORBH* 1 0 " 3 + $RE) "2+ ( $RE+$A236*1000) "2- (2*$ORBH*10"3 + 
2 * $RE) * ($RE + SA236*1000) * ©COS ( ©RADIANS ( $B2 3 6 - SNADIR) ) ) /1000 
CoherenceL : D2 3 6 : +SC235-SC236 

CoherenceL : E236 : © IF ( $A2 3 6>2 0 , 0 , © IF ( $A2 3 6<0 . 0 1 , 7 . 36*10" ( - 14 ) , 

©IF(SA236<3, 1 .59 *10" -12* ($A236*1000) " (-4/3) , 3 .66*10" (-17) ) ) ) 

CoherenceL : F2 36 : ©IF ( $C2 3 6>$R, 0 , $D2 3 6 * 10 " 3 *$E236 * ( 1-$C236/ $R) " (5/3) ) 
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Coher enceL : G2 3 6 : ( $ HCOHERENCE * (2*@PI/ ( $ LAMBDA* 10 A ( -6 ) ) ) n 2* 

@SUM ( F23 6 . . $F$237 } ) * (-3/5) 

CoherenceL : A2 3 7 : 20 

CoherenceL : B2 3 7 : ©DEGREES (©AS IN ( (ORBH*10~3+RE) * ©SIN (©RADIANS (NADIR) ) / 
(RE+$A237*1000) ) ) 

CoherenceL :C2 3 7 : @SQRT( ( $ORBH*10~3+$RE) A 2+ ( $RE+$A237*1000) ~2- (2*$ORBH*10 A 3+ 
2*$RE) * ($RE+$A237*1000) * ©COS ( ©RADIANS ( $B2 37- $NADIR) ) ) /1000 
CoherenceL : D237 : +$C236-$C237 

CoherenceL :E2 3 7 : ©IF ( $A237>20 , 0 , ©IF ( $A237<0 . 01 , 7 . 36 *10~ ( -14 ) , 

©IF ( $A237<3 , 1 . 59 * 10 ^ 12 * ($A237*1000) A (-4/3) , 3 . 66*10 A (-17) ) ) ) 

CoherenceL :F2 3 7 : ©IF ( $C237>$R, 0, $D2 37*10^3*$ E2 37* ( 1-$C237 / $R) A (5/3) ) 
CoherenceL :G2 3 7 : { $HCQHERENCE* (2*@PI/ ( $ LAMBDA* 10 ~ ( - 6 ) ) ) ~2* 

©SUM ( F237 ,.$F$237))~(-3/5) 


The following lines set up a search criteria table and determine the refractive structure function at 
the altitude of interest given on the Results page. 

CoherenceL :A23 9 : 'Search Criteria 

CoherenceL :A240 : 'Altitudel 

CoherenceL : C240 : 'Result 

CoherenceL :A2 41 : +ALTITUDE1>=$ALT/1000 

CoherenceL :C241 : ( SHCOHERENCE* (2*@PI/ ( $ LAMBDA* 10* ( - 6 ) ) ) / '2*@DSUM(WPL37- 

DAT, 5 , A240 . . A2 41) > ^ ( -3 / 5 ) 


4.9 Alignment 


The following code lists the contents of the Alignment page. The user sees a plot of loss due to 
misalignment. Columns A and B contain a parameterised misalignment loss curve. Column B 
contains the misalignment loss as a function of the parameter, aD (Jt/X) where a is the mis- 
alignment angle between the signal and BPLO, D is the optical beam diameter and X is the lidar 
operating wavelength. Cells D37 - E44 contain a search table to find the two values of aD (Jt/X) 
which straddle the misalignment angle given on the Instrument page. A linear interpolation 
between these two values is then performed to obtain the misalignment loss. Columns L and M 
contain the misalignment loss (column M) as a function of angle (column L) for the wavelength 
given on the Instrument page. These are the values plotted in the graph shown on this page. 


Alignment ;*A3 6 
Alignment : B3 6 
Alignment : L3 6 
Al ignment : M3 6 
A1 ignment : A3 7 
Al ignment : B37 
Al ignment : D37 
Al ignment : L37 
Al ignment : M37 
Al ignment : A3 8 
Alignment : B3 8 
Al ignment : D3 8 
Alignment : E3 8 
Alignment : G3 8 
Al ignment : L3 8 
Al ignment : M3 8 
Ai ignment : A3 9 
Al ignment : B3 9 
Al ignment : D3 9 
Al ignment : E3 9 
Alignment : G3 9 


' param 
' eta_align 
' del taangl 
' Loss 
0 
1 

'Search Criteria 
+ PARAM * $ LAMBDA / (!2PI*$D) 
1 - ETA_ALIGN 
0 . 1 

0 . 9987 

A Lower Limit 
A Upper limit 
'XY Data Pairs 
+ A 3 8 * $ LAMBDA / (@PI*$D) 

1 - 33 8 

0.2 

0 . 9951 
"param 
"param 
"Xlower 
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Alignment :H3 9 : '"Xupper 

Alignment : 139 : ^Ylower 
Alignment : J3 9 : * Yupper 
Alignment :L3 9: +A3 9 * $LAMBDA/ (@PI*$D) 

Alignment : M3 9 : 1-B39 
Alignment :A40 : 0.3 
Alignment : B40 : 0.9892 

Alignment :D40 : @ABS ( PARAM-DUMMY1 ) <0 . 1 # AND #PARAM<= DUMMY 1 

Alignment : E40 : +PARAM>DUMMY1 

Alignment :G40 : (3VLOOKUP ( DUMMY1 , ALIGNDAT , 0 ) 

Alignment :H40 : @DMIN (ALIGNDAT , 0, E39 . .E40) 

A1 ignment : 1 4 0 : @VLOOKUP ( DUMMY1 , ALIGNDAT # 1 ) 

Alignment : J40 : @DMAX ( ALIGNDAT , 1, E39 . . E40) 

Alignment : L40 : +A40 * $LAMBDA/ < @PI * $D) 

Alignment : M40 : 1-B40 

Alignment : A41 : 0.4 

Alignment : B41 : 0.9812 

Alignment : L41 : +A41 *$ LAMBDA/ (@PI*$D) 

Alignment :M41 : 1-B41 

Alignment : A4 2 : 0.5 

Alignment : B4 2 : 0.9708 

Alignment :L42 : +A4 2 * $LAMBDA/ (@PI*$D) 

Alignment : M4 2 : 1-B42 

Al ignment : A4 3 : 0.6 

Alignment : B4 3 : 0.9583 

Alignment : D4 3 : "Search value (X) = 

Alignment : E43 : @ABS ( ALIGN_ANGLE ) * D * @ P I / LAMBDA 
Alignment :L43 : +A43 * $ LAMBDA/ (@PI*$D) 

Alignment : M4 3 : 1-B43 
Alignment : A4 4 : 0.7 
Alignment : B4 4 : 0.9436 
Alignment : D44 : "Solution (Y) = 

Alignment : E44 : @IF (DUMMY1 = 0 , 1 , @IF {G40=H40 , 140, +140+ ( J40-I40) / (H40- 
G40 ) # ( DUMMY1-G40 ) ) ) 

Alignment :L44 : +A44 * $LAMBDA/ (@PI*$D) 

Al ignment : M44 : 1-B44 

Alignment : A4 5 : 0.8 
Alignment : B4 5 : 0.927 
Alignment : L45 : +A45 * $LAMBDA/ (@PI*$D) 

Alignment : M4 5 : 1-B45 
Al ignment : A4 6 : 0.9 
Alignment : B4 6 : 0.9084 

Al ignment : L4 6 : +A4 6 * $ LAMBDA/ ( @PI * $D) 

Alignment : M4 6 : 1-B46 

Al ignment : A47 : 1 

Al ignment : B47 : 0.8881 

Al ignment : D47 : 'This page searches the misalignment database to find the val- 

ues that straddle 

Alignment :L47 : +A47 * $LAMBDA/ (@PI*$D) 

Al ignment : M47 : 1-B47 

Al ignment : A48 : 1.1 
Al ignment : B48 : 0.8661 

Alignment : D48 : 'the angular misalignment. A linear interpolation is then used 

between these two 

Alignment :L4 8 : +A48 # $ LAMBDA/ (@PI*$D) 

Alignment : M4 8 : 1-B48 

Al ignment : A49 : 1.2 

Al ignment : B49 : 0.8427 

Al ignment : D49 : 'points to determine the misalignment loss. The use of a linear 

interpolation 

Alignment : L4 9 : +A49 * $LAMBDA/ <3PI*$D) 

Al ignment : M4 9 : 1-B49 

Al ignment : A50 : 1.3 

Al ignment : 850 : 0.8179 
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Alignment : D50 : 
imal place this 
Alignment : L50 : 
Alignment : M50 : 
Alignment : A51 : 
Alignment : B51 : 
Alignment : D51 : 
Alignment : L51 : 
Alignment : M51 : 


'introduces an error in the misalignment loss in the third dec- 
is 

+A50 * $LAMBDA/ (@PI*$D) 

1-B50 

1.4 

0.7919 

' insignificant . 

+A5 1 * $LAMBDA/ (@PI*$D) 

1-B51 


The formulae in columns A,B.L and M are repeated until the data for the plot ( columns L and M) 
are completed:- 


Alignment : A156 
Alignment : B156 
Alignment : LI 5 6 
Alignment : Ml 5 6 
Alignment : A157 
Alignment : B157 
Alignment : LI 5 7 
Alignment : Ml 5 7 


19 . 5 

0 .00023224357861854 
+A156*$LAMBDA/ (@PI*$D) 
1-B156 
20 

0.0002147988979521 
+A1 57*$ LAMBDA/ (@PI*$D) 
1-B157 


The remaining rows then contain the rest of the misalignment look-up table until:- 

Alignment : A172 : 95 

Alignment :B172: 1 . 9573 1733773E-06 

Alignment :A17 3 : 100 

Alignment : B17 3 : 1 . 67746225269E-06 


4.10 Performance 

The Perf page contains two plots, one of the wideband and narrowband SNR as a function of 
backscatter and one of the velocity error and probability of correct detection as a function of back- 
scatter. The SNR plot contains a button which causes these graphs to be updated. The data for 
these plots is generated by the macro given in and is placed in cells A36 E59. Column A contains 
the backscatter value, B/C the wideband and narrowband SNR respectively and D and E contain 
the probability of a good estimate and the velocity error respectively. 

Unlike the other plots which update themselves automatically when one of the relevant lidar 
parameters changes, these plots only update when the user clicks on the “Update Graphs’* button. 
This will be discussed later. 

This is the last page of direct relevance to the casual user - the remaining pages contain coding 
and although currently visible to the user they could be made invisible to prevent accidental alter- 
ation. 


4.11 Constants 

The following code lists the contents of the cells on the Constants page:- 
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Constants 

A2 

'h 


Constants 

B2 

662.60755 

*10*- 

Constants 

C2 

' J/s 


Constants 

A3 

'c 


Constants 

B3 

299792458 


Constants 

C3 

' m/s 


Constants 

A4 

're 


Constants 

B4 

6371315 


Constants 

C4 

' m 


Constants 

A5 

' me 


Constants 

B5 

5.979*10* 

24 

Constants 

C5 

'kg 


Constants 

A6 

' te 


Constants 

B6 

8.639988* 

10*4 

Constants 

C6 

' sec 


Constants 

A7 

'G 


Constants 

B7 

6.67259*10^-11 

Constants 

A8 

1 hcoherence 

Constants 

B8 

2.914383 


Constants 

A9 

"1/pulse 

lngth 

Constants 

B9 

0 . 187 4/TAU 

Constants 

C9 

'MHz 



4.12 Maximum liklihood estimator 

This page contains a look-up table for the maximum liklihood estimator. The table is shown 
below together with the column and row assignments. 



B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

l 

Omega 

M 

bO 

alpha 


chi 

go 

epsilon 

delta 

mu 

2 

3 

20 

4 . 6308 

1 . 1923 

1.2535 

0.9145 

6.2318 

1.3798 

0.5307 

0.1462 

3 

3 

25 

6 . 6127 

1 . 1843 

1.573 

0.887 

4.3866 

2.2939 

0.2345 

0.1082 

4 

3 

30 

9 . 6148 

1.18 

2 . 0453 

0.9001 

3.5357 

3.1954 

0.1401 

0.0807 


3 


13 .067 

1 . 1495 

2 . 5687 

0.8758 

4.2848 

2.5322 

0.2026 

0.1 


3 

40 

14 . 003 

1 .1862 

2.7099 

0.8854 

3.5066 

4.0879 

0.1054 

0.0761 

7 

6.7 


0.2318 

4.0951 

0.0196 

0.2138 

17.454 

1 . 1469 

0.8738 

0.4349 

8 

6.7 


0.7783 

5 .6302 




0.9637 

1.2103 

0.3021 

9 

6 . 7 


2 . 2702 

1.5914 

0.4216 

0.7402 

12.063 

1.1882 

0.8432 

0.2037 

10 

6 . 7 


5.7615 

1 . 1484 

1.1356 

0 . 8928 

9 .3993 

1.2744 

0.708 

0.1468 

11 

6 . 7 

40 

9.832 

1.0981 

1 .8154 

0.9536 

7.7512 

1.4531 

0.544 

0.1122 | 


4.13 Lists 

This page contains small look-up tables used by the program. The tables are shown below, 
together with their row/column assignments. The first table is a list of wavelengths for which 
there are entries in the extinction database spreadsheet, ATMOS.WB1. This table appears as the 
pick list in the laser wavelength entry box. 
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A 

1 

2.012552 

2 

2.021842 

3 

2.065479 

4 

2.091282 

5 

9.1145187 

6 

9.21969 

7 

9.249946 

8 

10.591043 


The following table lists the available receiver geometries and their beam truncation factors. 



C 

D 

E 

F 

G 

i 

Geometry 

Tt 

etah 

gammar 


2 

Wang 

0.955 

0.42 

0.802 

0.802 

3 

Matched Rye 

0.929 


0.87 

0.87 

4 

Unmatched Rye 

0.951 

0.461 

1.186 

0.815 


This small table is used in the receiver selection box to choose between complex and real receiver 
geometries (this affects the signal processing). 



I 

1 

Complex 

2 

Real 


The next table lists the available scanning options and is used by the scanner input dialog. 
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J 

1 

Dual-look 

2 

Conical 

3 

Wedge 


4.14 Macros 

The following table is a coded macro to calculate and plot the performance plots (Section 4. 10). It 
is executed when the user clicks on the “Update Graphs” button on the SNR plot. Column A con- 
tains variable names and labels, column B contains the actual coding or the variable contents and 
column J contains comments on the macro steps. 



A 

B 

J 

2 

THIS MACRO RECALCULATES THE INSTRUMENT PERFORMANCE AS A FUNCTION OF 
BACKSCATTER 

3 

AND UPDATES THE PERFORMANCE GRAPHS SO THAT ONLY PERFORMANCE FOR VAL- 
UES OF 

4 

PROBABILITY > mm_prob ARE PLOTTED. THIS IS BECAUSE THE VELOCITY ESTIMATION 
THEORY IS 

5 

UNCERTAIN FOR PROBABILITIES LESS THAN THIS. IN PRACTICE IT IS UNLIKELY A 
LIDAR OF SUCH POOR 

6 

SINGLE SHOT PERFORMANCE WOULD BE BUILT! 

7 




8 

save_beta 

IE-06 


9 

beta_count 

24 


10 

mm_beta_point 

50 


1 1 

min_prob 

0.2 


12 




13 

perform _c ale 

{ EDITGOTO perform: A 1 1 1 LET $save_beta. 
SB ETA 1 

Store beta value 

14 


\ FOR beta_count,0,23, 1 ,beta_calc ) 

Repetitive calculation loop, 
increments beta value pointer 

15 


| LET $BETA,$save_beta) 

Restore beta to original value 
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{WINDOWSOFF}) GraphEdit SNR} 

Turn window updating off and 
edit the SNR graph 

{ SERIES. Data_Range L+" Perform :B"&@- 
STRING(min_beta_point,0)&"..B59"} 

Change graph series to match 
new performance 

{ SERIES. Data_Range 2,VTerform:C ,, &@- 
STRING(min_beta_point,0)&"..C59"} 

Only performance for values 
of probability > 0.2 are plotted 


( SERIES. Data_Range , XAxisLabelSeries”,+ M - 
Perform:A'&@STRING(min_beta_point 1 0)&". 
A59" } 


{ Series. Go } { WindowClose } 


(GraphEdit Vel_error) 


( SERIES. Data_Range l,+"Perform:- 
E’&@STRING(min_beta_point,0)& M ..E59 M } 


( SERIES. Data_Range 2,+ 'Perform:D ,, &@- 
STRING(min_beta_point,0)&"..D59 M } 


{ SERIES. Data_Range , XAxisLabelSeries M ,+ M - 
Perform:A , ’&@STRING(min_beta_point,0)& M 
A59 M } 


{ Series. Go } { WindowClose } 



Instigate graph changes and 
close edit window 


Edit velocity error/probability 
graph 


Change graph series to match 
new performance 


Only performance values for 
probability > 0.2 are plotted 



Instigate graph changes and 
close edit window 


27 

beta_calc 

( Let SB ETA, @INDEX( performance ,0,beta_- 
count)( 

Set beta to value in perfor- 
mance table 

28 


— — 

(CALC) 

Calculate performance 

29 


(Let + perform:B M &@STRING(beta_- 
count+36,0),snr_wide ) 

Store wideband SNR in per- 
formance table 

30 


(LET +" perform :C'&@STRlNG(beta_- 
count+36,0),snr_narrow ) 

Store narrow band SNR in 
performance table 

31 


(LET + "perform :D”&@ STRING(beta_- 
count+36,0).pgood ) 

Store probability of good esti- 
mate in performance table 

32 


( LET +’ "perform: E”&@STR ING(beta_coun- 
t+36,0),sigmav_i \ 

Store velocity error in perfor- 
mance table 


33 


1 IF pgood<mm_prob 1 { LET min_beta _poinL(be- 
ta_count+36)( 


If probability is < 0.2 incre- 
ment minimum beta pointer. 



































4.15 Other pages 


The spreadsheet also contains two other pages - Messages and Block_Names. The Messages page 
will contain help messages which are accessible to the user (not currently implimented) whilst the 
Block_Names page contains a tabulation of the variable names assigned to cell locations. This 
page is not required for operation of the program but is included to assist in the further develop- 
ment of the program. 
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5.0 Appendix D - ANSI Z136. 1-1992 Eye Safety Calculation 

5.1 Program structure 

Program EYES92 - main routine, allows user to select the MPE (Maximum Permissible Expo- 
sure) as 

• a function of wavelength for a fixed pulse length or 

• a function of pulse length for a fixed wavelength or 

• a function of velocity resolution over a wavelength range (approximation). 

The program saves the results to the data file mpevsl.dat (mpe vs wavelength) or mpevst.dat (mpe 
vs. pulse length). The output is in two columns, the first contains the variable (pulse length (s) or 
wavelength (pm) and the second contains the MPE (J/cm 2 ). 

Function MPE92(LB,T) - is the function that calculates the MPE for a given wavelength, LB in 
microns and a given pulse length, T in seconds. Called by the main program. 

Subroutine CLRSCRN(I) - is a utility subroutine that scrolls the output screen by I lines. Called 
by the main program. 


5.2 Program code 
5.2.1 Program EYES92 


c 


c 


PROGRAM EYES92 


IMPLICIT NONE 

REAL* 4 RESL , MIN , MAX , STEP , VAL , RESULT , LB , TP , MPE92 , T 
INTEGER* 2 OPTION 
CHARACTER* 11 FILNAM 


CALL CLRSCRN 
WRITE ( * , * ) ' 

WRITE ( * , *) ' 

WRITE ( * , * ) ' 

WRITE {*,*) ' 

WRITE ( * , * ) ' 

WRITE ( * , * ) ' 

WRITE { * , * ) ' 

WRITE { * , *) 
WRITE (*,*) ' 

WRITE (*,*) ' 

WRITE { * , * ) ' 

WRITE { * , * } ' 

WRITE ( * , * ) ' 

WRITE { * , * ) ' 

WRITE ( * , * ) ' 

WRITE ( * , * ) ' 

WRITE ( * , * } ' 

WRITE ( * , * ) 
PAUSE 'Type 


(40) 

1992 ANSI STANDARD EYE SAFETY CALCULATIONS ' 
Gary D. Spiers' 

Center for Applied Optics' 

University of Alabama in Huntsville' 
Huntsville ' 

AL 35899' 

(205) 895 6030 ext. 448' 

This program calculates the maximum permissible' 
exposure (MPE) in J/crrT2 for direct ocular exp-' 
osure, intrabeam viewing to a laser beam. The ' 
valid wavelength range is 0.4-1000 microns and' 
the valid pulse durations (half amplitude po-' 
ints) are Ins to 3000s. Pulses shorter than Ins' 
are treated as Ins pulses at present due to 1 
inadequate safety data. All calculations are ' 
based on the ANSI Standard Z136 . 1- 1992 . ' 

the return key to continue' 
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CALL CLRSCRN { 40 ) 

100 WRITE ( * , * ) 

WRITE (*,*) 'Three options are available:-' 

WRITE ( * , *) 

WRITE (*,*) '1) MPE as a function of wavelength for a ' 

WRITE (*,*) 1 fixed pulse length. ' 

WRITE ( * , * ) ' OR' 

WRITE (*,* ) '2) MPE as a function of pulse length for a ' 

WRITE (*,* ) ' fixed wavelength.' 

WRITE (*,*) ' OR' 

WRITE ( * , * ) '3) MPE as a function of velocity resolution' 

WRITE {*,*) ' over a wavelength range.' 

WRITE ( * , *) 

WRITE (*,*) 'ENTER OPTION NUMBER:' 

READ ( * , * ) OPTION 
WRITE ( * , *) 

IF ((OPTION . NE . 1 ) .AND. { OPTION . NE . 2 ) .AND. ( OPTION . NE . 3 ) ) THEN 
WRITE {*,* ) 'Please use a valid option number!!!!' 

GOTO 100 

ELSEIF (OPTION. EQ . 1 ) THEN 
WRITE ( * , * ) 'ENTER PULSE LENGTH (MICROSECONDS):-' 

READ ( * , * ) TP 
T=TP* 0 . 000001 
IF ( T .LT. IE-9) THEN 
WRITE ( * , *) 

WRITE (*,*) 'Pulselength set to Ins.' 

T = IE-9 
END IF 

IF (T .GT. 3000) THEN 
WRITE {*,*) 

WRITE (*,*) 'Pulse length set to 3000s.' 

T =3000 
END IF 

ELSEIF (OPTION . EQ . 2) THEN 
130 WRITE ( * , * ) 'ENTER WAVELENGTH (MICRONS) 

READ ( * , * ) LB 

IF ((LB .LT. 0.4) .OR. (LB .GT. 1000)) THEN 
WRITE ( * , *) 

WRITE ( * , * ) ‘The wavelength must be in the range 
@ 0.4-1000 microns.' 

GOTO 130 
END IF 

140 WRITE (*,*) 'ENTER THE MINIMUM AND MAXIMUM PULSE' 

WRITE (*,*) 'LENGTHS OF INTEREST (microseconds) 

READ ( * , *) MIN, MAX 

IF ( (MIN. LT. 0.001) .OR. ( MAX . GT . 3 0 00 ) .OR. (MAX . LT . MIN) ) THEN 
WRITE(*,*) ' Please provide valid values in the range' 

WRITE ( * , * ) 1 0.001 to 3000 microseconds!' 

GOTO 140 
END IF 

ELSEIF (OPTION . EQ . 3) THEN 
150 WRITE ( # ,M 'ENTER VELOCITY RESOLUTION (m/s) 

READ ( * , * ) RESL 
IF (RESL .LE. 0.0) THEN 
WRITE ( * , # ) 

WRITE { * , * ) ‘The velocity resolution must be a 
@ positive number.' 

GOTO 150 
ENDIF 
END IF 

IF (OPTION .NE. 2) THEN 
WRITE ( * , * ) 

WR I T E ( * , * ) 'Now set up wavelength range of interest.' 

WRITE (*,*) 'Wavelength increment is 5nm. ' 

200 WRITE ( * , * ) 
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WRITE ( * , * ) ' ENTER MINIMUM AND MAXIMUM WAVELENGTH ' 

WRITE ( * , *) ' (MICRONS) 

READ ( * , * ) MIN, MAX 
WRITE ( * , *) 

IF ( (MIN . LT .0.4) .OR. (MAX . GT . 1000 ) .OR. (MAX . LT . MIN) ) THEN 
WRITE(*,*) ' Please provide valid values in the range' 
WRITE(*,*) ' 0.4 to 1000 microns !!!!! ' 

GOTO 200 
ENDIF 

WRITE (*,* ) 'The data will be written to the file mpevsl.dat.' 
F I LNAM= ' mpe vs 1 . da t ' 

ELSEIF (OPTION . EQ . 2) THEN 

WRITE { * , * ) 'The data will be written to the file mpevst.dat.' 
WRITE (*,* ) 'The pulse length increment is 5 ns.' 

FILNAM= ' mpevst . dat ' 

ENDIF 

OPEN ( 1 , FILE=FILNAM) 

STEP=0 . 005 

DO 300 VAL = MIN, MAX, STEP 
IF (OPTION .EQ. 1) THEN 
LB=VAL 

ELSEIF (OPTION . EQ . 2) THEN 
T=VAL* 0 .0 00001 
ELSEIF (OPTION . EQ . 3) THEN 
LB = VAL 

T=LB* IE- 6 / ( 2 *RESL ) 

ENDIF 

RESULT=MPE92 (LB, T) 

IF (OPTION . NE. 2) THEN 
WRITE (1, *) LB, RESULT 
ELSEIF (OPTION . EQ . 2) THEN 
WRITE (1,*) T, RESULT 
ENDIF 

300 CONTINUE 
STOP 
END 


5.2.2 Function MPE92 


REAL *4 FUNCTION MPE92(LB,T) 

REAL*4 CA, CB, T1 , LB, CC, T 

C 

IF ((LB .GE. 0.4) .AND. (LB . LT . 0.7)) THEN 
CA= 1 

ELSEIF ((LB .GE. 0.7) .AND. (LB . LT . 1.050)) THEN 
CA=10** (2 . OMLB-O .7) ) 

ELSEIF ((LB .GE. 1.050) .AND. (LB . LT . 1.4)) THEN 
CA=5 
ENDIF 

C 

IF ((LB .GE. 0.4) .AND. (LB . LT . 0.55)) THEN 
CB = 1 

ELSEIF ((LB .GE. 0.55) .AND. (LB . LE . 0.7)) THEN 
CB=10* * (15* (LB-0 . 55) ) 

ENDIF 

C 

IF ((LB .GT. 1.05) .AND. (LB . LE . 1.15)) THEN 

CC=1 . 0 

ELSEIF ((LB .GT. 1.15) .AND. (LB . LE . 1.2)) THEN 
CC=10»* (18* (LB-1 . 15) ) 

ELSEIF ((LB .GT. 1.2) .AND. (LB . LE . 1.4)) THEN 
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c 


CC=8 

ENDIF 

IF ((LB .GE. 0.55) .AND. (LB .LE. 0.7)) THEN 
Tl=10*10** (20* ( LB- 0 . 55) ) 

ENDIF 

IF (LB .LE. 0.7) THEN 
IF (T .LE. 1.8E-5) THEN 
MPE92=5E-7 

ELSEIF ( ( T .GT. 1.8E-5) .AND. (T . LT . 10)) THEN 
MPE92=1 . 8*T**0 . 75*lE-3 


ELSEIF ((LB .LT. 0.55) .AND. 
@ . LE . 1E4 ) ) THEN 

(T 

.GT. 

10) 

.AND. 

(T 



MPE92=0 .01 

ELSEIF ((LB .GT. 0.55) .AND. 

(T 

.GT. 

10) 

.AND. 

(T 

.LE. 

Tl) 

<3 ) THEN 








MPE92=1 . 8*T**0 . 75*lE-3 
ELSEIF ((LB .GT. 0.55) .AND. 

(T 

.GT. 

Tl) 

.AND. 

(T 

.LE. 

IE 

@4)) THEN 








MPE92=10*CB*0 .001 
ELSEIF ((LB .GE. 0.4) .AND. 

<T . 

.GT. 

1E4 ) 

.AND. 

<T 

.LE. 

3E 


@4)) THEN 

MPE92=CB*lE-6 

ENDIF 

ELSEIF ((LB .GT. 0.7) .AND. (LB . LE . 1.05)) THEN 
IF (T .LE. 1.8* IE-5 ) THEN 
MPE92=5*CA*lE-7 

ELSEIF ( (T .GT. 1.8*lE-5) .AND. (T .LE. 1E3)) THEN 
MPE92=1 . 8*CA*T**0 .75*0 .001 
ELSEIF ( (T .GT. 1E3) .AND. (T .LE. 3E4)) THEN 
MPE92=3 20*CA* IE- 6 
ENDIF 

ELSEIF ((LB .GT. 1.05) .AND. (LB .LE. 1.4)) THEN 
IF (T .LE. 5E-5) THEN 
MPE92=5*CC*lE-6 

ELSEIF ( (T .GT. 5E-5) .AND. (T .LE. 1E3)) THEN 

MPE92=9*CC*T**0 . 75*0. 001 
ELSEIF ( (T .GT. 1E3) .AND. (T .LE. 3E4)) THEN 
MPE92= 1 . 6 *CC* IE- 3 
ENDIF 

ELSEIF ((LB .GT. 1.4) .AND. (LB . LE . 1.5)) THEN 

IF (T .LE. IE-3 ) THEN 
MPE92=0 . 1 

ELSEIF ((T .GT. IE-3) .AND. (T . LE . 10)) THEN 

MPE92=0 . 56*T**0 . 25 

ELSEIF ( ( T .GT. 10) .AND. (T . LE . 3E4)) THEN 

MPE92=0 . 1 
ENDIF 

ELSEIF ((LB. GT. 1.5) .AND. (LB .LE. 1.8)) THEN 

IF (T .LE. 10) THEN 
MPE92=1 

ELSEIF ( (T .GT. 10) .AND. (T .LE. 3E4)) THEN 
MPE92=0 . 1 
ENDIF 

ELSEIF ((LB .GT. 1.8) .AND. (LB . LE . 2.6)) THEN 
IF (T .LE. IE- 3 ) THEN 
MPE92=0 . 1 

ELSEIF (( T .GT. IE-3) .AND. (T . LE . 10)) THEN 
MPE92=0.56*T**0.25 

ELSEIF ( (T .GT. 10) .AND. (T .LE. 3E4)) THEN 
MPE92=0 . 1 
ENDIF 

ELSEIF ((LB .GT. 2.6) .AND. (LB . LE . 1E3 ) ) THEN 

IF (T .LE. IE-7) THEN 
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(T .LE. 10)) THEN 


MPE92=lE-2 

ELSEIF {(T .GT. IE-7) .AND. 

MPE92=0.56*T**0.25 
ELSEIF ( (T .GT. 10) .AND. (T . LE . 3E4)) THEN 
MPE92=0 . 1 
ENDIF 
ENDIF 
RETURN 
END 


5.2.3 Subroutine CLRSCRN 


SUBROUTINE CLRSCRN ( I ) 
INTEGER J, I 
DO 100 J=1 ,1,1 
WRITE ( * , *) 

100 CONTINUE 
RETURN 
END 


5.3 Sample plot of data 


Eye Safety Limits 
(0.5 ps pulse length) 
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6.0 Appendix E - Eye Safety Analysis 

The following paper was presented as Paper TuA6 at the 7th. Coherent Laser Radar Applications 
and Technology Conference, Paris, France, 19th-23rd July 1993. 
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Eye Safety Considerations For Selecting The Wavelength of a Space Based 

Coherent Doppler Wind Lidar 
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Center for Applied Optics 
University of Alabama in Huntsville 
Huntsville, AL 35899, USA 
(205) 895 6030 ext. 448 

Michael J. Kavaya 
NASA Marshall Space Flight Center 
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David A. Bowdle 
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Introduction 

This paper derives an expression for the eye safety of a coherent doppler wind lidar starting 
from the signal to noise ratio (SNR) equation. The intensity at the eye of a potential observer and the 
proposed 1993 ANSI eye safety standard are folded in to produce an eye danger index. 


Theory 

The SNR from an aerosol target (subscript a) at range, R a , from a lidar operating at a 
wavelength, X, is given by: 


SNR (A.) = 


nD^U L Tl (X) K u (X) K „ (X) cp a (X) T, (X) T r (X) 


8hl)B (X) R, 


1 + 


' °b y 
2P 0 .(*)J. 


( 1 ) 


where U^ is the laser pulse energy at the laser, K^X) and K^X) are the atmospheric transmission of 
the transmit and receive beams, c is the speed of light, 0 a (X) is the atmospheric backscatter 
coefficient, T t (X) and T r (X) are the transmittance of the transmit and receive optics, h is Planck’s 
constant, v is the laser frequency, B(X) is the processing noise bandwidth and Po a (X) is the transverse 
field coherence length. f|(X) is the efficiency of the lidar which is given by: 

= n s (^)Ti a (X)ti p (X)n q (>»)Ti u (X) (2) 

where T1 S (X), q a (X), rj p (X), q q (X) and T1 U (X) are the lidar system, misalignment, polarisation, 
effective heterodyne quantum and unexplained loss efficiencies respectively. For a Gaussian beam 
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the 1/e 2 intensity diameter of the transmitted beam, Dg, is optimised relative to the telescope 
physical diameter, Dj, when [1]: 

D b = 0.802 D t (3) 

The processing noise bandwidth, B(X) is: 

2V s 

B (X) = (4) 

where V s is the velocity search window. 

The intensity, Qg at the eye of an observer (subscript e) at a distance, R^, from the lidar, 
including diffraction but assuming negligible spreading of the beam by refractive turbulence is given 
by: 


Q e W = 


4nD B U L K te (X) T, (X) M 
16R 2 X 2 + n 2 Dg 


(5) 


where K te (X) is the atmospheric transmission to the observer and M is a factor to account for the 
transverse intensity distribution of the optical beam. For example, for a fixed pulse energy, the 
intensity on the axis of a Gaussian beam is twice that of a uniform beam with the same diameter and 
therefore M=2 for a Gaussian profile. 

An eye danger index X(X) can be defined by: 


X(X) 


Q e W 

QmpeW 


( 6 ) 


where Q mpe (X) is the proposed 1993 ANSI eye safety standard maximum permissible exposure 
(MPE) value shown in Figure (la) for a lps pulse length. For values of X(X) > 1 the lidar exceeds 
the ANSI specification. Combining Equations (1,2, 4-6) gives: 


X(X) 


16SNR(X)hV s MR 2 K te (4pQ t (X) +D%) 

n (X) K rl (X) K u (X) T r (X) X 2 P t (X) Q mpc (X) p 0t (X) 2 (n 2 Dj + 16X 2 R 2 ) 


Figures (lb) and (lc) show the value of X(X) over the wavelength range 1-12 pm for a space 
based lidar and an earth based observer using the parameters listed in Table (1). The calculations 

assume a spherical earth, a lidar altitude of 525 km, a 45° nadir scan angle, an aerosol target at an 
altitude of 10 km and an observer on the ground.. The minimum detectable SNR was taken as -13.2 

dB, independent of wavelength. Both a X’ 2 (Figure (lb)) and a X'* (Figure (lc)) dependence of 
backscatter on wavelength were used to represent typical values near 10 km altitude. The 

backscatter was anchored to a value of 5 x 10* 11 (m-sr)* 1 at 9.11 |im, again a typical value for 10 
km. A laser pulse length of 1 ps was assumed when determining the MPE values. For wavelengths 
shorter than 2.6 pm the MPE does not depend on pulse length for pulses of width 1 ns - 50 ps, and 
only a weak pulse length dependence exists for wavelengths longer than 2.6 pm. The transverse 

coherence length has a X^ dependence and was anchored to a value of 58 m at 9.11 pm. The 
wavelength dependent atmospheric transmission is not suited to a simple expression and a constant 
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value of 0.86 was assumed for the atmospheric transmission to the observer. Atmospheric extinction 
to the aerosol target was assumed to be negligible. 


SNR(X) 

-13.2 dB 

h 

6.62607 xlO' 34 Js 

v s 

50 ms' 1 or ± 25 ms -1 

Tr(X) 

0.95 

M 

2 

W 

5 x 10' n /(m-sr) <ffi 9.11 nm 

Qmpe(^) 

ANSI standard (Figure la) 

R a 

754.2 km 

Re 

775.9 km 

Dj 

1.5 m 


w 

0.46 

n<A) 

0.5 

Tlp(A) 

1 

Tlq(X) 

0.4 

T1„(X) 

0.5 

K ra (X) 

1 

K to (X) 

1 

K te (X) 

0.86 

P0a(X) 

58m@9.U|im 


Table 1: Parameter values used to calculate Figures (lb) and (lc). 


Discussion 

We have combined the coherent lidar SNR equation with the ANSI eye safety standard to 
arrive at an expression for eye danger index. Equation (7) is valid for any location of a coherent 
Doppler lidar, aerosol target and human observer; however Figures (lb) and (lc) apply only to the 
specific space based configuration described in Table 1. For the space based configuration, the 
parameters in Equation (7) with the largest dynamic range are X, Q mpc (X) and (3(X). Limitations of 
the existing calculations are the use of wavelength independent terms for misalignment loss, 
atmospheric transmission, SNR threshold, detector quantum efficiency and exponent for the 
backscatter wavelength dependence. 
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Wavelength, (jim) 


Figure la) The MPE for a 1 ns pulse, b) and c) The eye danger index for the two backscatter 
dependencies on wavelength discussed in the text. 
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